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Abstract 

We present the analytic calculation of the McUin moments of the structure functions i^2, ^3 
and Fl in perturbative QCD up to second order corrections and in leading twist approximation. 
We calculate the 2-loop contributions to the anomalous dimensions of the singlet and non-singlet 
operator matrix elements and the 2-loop coefficient functions of F2, ^3 and Fl- We perform the 
inverse Mellin transformation analytically and find our results in agreement with earlier calculations 
in the literature by Zijlstra and van Ncerven. 



1 Introduction 



As one of the best studied reactions today, deep inelastic lepton-hadron scattering establishes the 
scale evolution of the structure functions, one of the most important precision tests of perturbative 
QCD. It provides unique information about the deep structure of the hadrons and most importantly, 
measurements of these structure functions in deep inelastic scattering (DIS) allow to extract the parton 
densities, which are subsequently used as input for many other hard scattering processes. 

Over the years, the ever increasing accuracy of deep inelastic and other hard scattering experiments 
has created a steady demand for more accurate theoretical predictions. The qualitative quest for 
understanding scaling violations of the structure functions in terms of asymptotic freedom in QCD [||] 
soon changed to the quantitative task of reducing theoretical uncertainties due to higher order QCD 
corrections and the necessary calculations to obtain next-to-leading order (NLO) perturbative QCD 
predictions were performed in refs.g |, |, |, |, |, |, |, 0, |ll|, ||]. 

Today, high precision analyses of experimental data, such as the determination of the strong 
coupling constant as and the parton densities call for complete next-to-next-to-leading order (NNLO) 
perturbative QCD predictions. This has motivated the calculation of the 2-loop coefficient functions 



of all DIS structure functions by Zijlstra and van Neerven 14, ^] and they could check their 
results for a number of fixed Mellin moments However, to complete the NNLO analysis of DIS 
structure functions one still has to know the perturbative QCD predictions at 3-loops for the anomalous 
dimensions of these structure functions and also the coefficient functions for the longitudinal structure 
function entering in the ratio of the longitudinal over the transverse cross section, R = (Jl/ctt- These 
quantities are still unknown, except for a small number of fixed Mellin moments [l^, 19 1, some of 
them related to sum rules. The results of refs.|19], have already been used in NNLO analyses ||l9| , pc| ] 
Unfortunately, this limited information about some fixed Mellin moments is generally not sufficient 
to allow for NNLO analyses of all data of DIS and related hard scattering experiments, such as the 
Drell-Yan process. As a consequence, there are still considerable uncertainties on the parton densities, 
a prominent example being the gluon density at small x. 

The aim of this paper is first of all to provide an independent check on the results of refs.|13, 14, |l^ ], 
by means of a completely different method. At the same time, we also wish to demonstrate the power 
of our method, which we believe to be flexible enough and most promising in view of the ultimate 
challenge, the calculation of the anomalous dimensions at 3-loops. Our approach, that actually dates 
back to the origins of QCD |l], ^ is to calculate the Mellin moments of the DIS structure functions 
analytically as a general function of A^. This idea was further pioneered by Kazakov and Kotikov |^] 
to obtain the longitudinal structure function Fi at two loops. In the present paper, we calculate in 
this way the Mellin moments of all unpolarized DIS structure functions F2, -F3 and up to two loops. 
Subsequently, we perform the inverse Mellin transformation to express our results in momentum space 
as functions of x to compare with refs.|p^, U, p^ . 

A remark about other strategies towards completing the NNLO perturbative QCD predictions for 
DIS structure functions is in order here. A straightforward extension of the methods of ref.||l^ to 
simply calculate sufficiently more fixed Mellin moments for a combined NNLO analyses of all data 
of, say, DIS and Drell-Yan experiments is not feasible. One has to know the Mellin moments at least 
up to = 100 which would result in analytical expressions consisting of huge rational numbers. 
Unfortunately, this is beyond the capabilities of present computer algebra programs. A different 
approach is the direct calculation of the anomalous dimensions of DIS operator matrix elements, 
extending the work of ref.[^, Q to 3-loops. First steps in this direction have been achieved in ref . , 
where the finite terms of DIS operator matrix elements at two loops have been calculated. 

The outline of this paper is as follows. In section |2| we briefiy discuss the operator product expansion 
(OPE) and some issues of the renormalization procedure. Section ^ gives a detailed explanation of the 
method to calculate the Mellin moments as an analytical function of A^ for a given diagram. It also 
contains a short summary of properties of harmonic sums |Q, 22, 23] and harmonic polylogarithms [p^ ]. 
Section ^ describes our calculation and lists our results for the structure functions F2, F3 and F^ up 
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to two loops, both in Mellin space as a function of N and in momentum space as a function of 
X expressed in terms of harmonic polylogarithms. Finally, section ^ gives our conclusions. The 
appendices contain some relations between the standard polylogarithms and Nielsen functions and 
the harmonic polylogarithms and all explicit expressions for the 2-loop coefficient functions. 



2 Formalism 

In this section we set up the stage for the calculation of Mellin moments of the deep inelastic structure 
functions. We briefly recall the operator product expansion and, in particular, pay attention to the 



renormalization procedure. For reviews see refs.|2^, 26|. 

We wish to calculate the hadronic part of the amplitude for unpolarized deep inelastic lepton- 
nucleon scattering which is given by the hadronic tensor 

W,u{p,q) = ^/d^ze'^-(P|jt(z)J.(0)|P) (1) 

1 1 p'^q^ 
= e^y —Fl{x, Q^) + d^y —F2{x, Q"^) + le^yap F^{x, Q^) , 

where is either an electromagnetic or a weak current and |P) is the unpolarized hadronic state. 
The boson transfers momentum q, = —q^, the hadron carries momentum p and the Bjorken scaling 
variable is defined as x = Q'^/ {2p ■ q) with < x < 1. The tensors e^jy and d^jy are given by 

1x 

df,u = -g,iu - PiJiPv-^ - {Pf^Qu + Pi'Qf^)-^ ■ (3) 
The longitudinal structure function Fl is related to the structure function Fi, 

Fl{x, Q2) = F2(x, Q2) _ 2xFi{x, Q^) . (4) 

The structure function F^ describes parity-violating effects that arise from vector and axial-vector 
interference. It vanishes for pure electromagnetic interactions. 

We are interested in the Mellin moments of the structure functions, defined as 

1 

Fi^iQ^) = J dxx^-^Fi{x,Q^), i = 2,L. (5) 



A similar relation defines F^ with replaced by -|- 1 in the integral on the right hand side in eq. (^) . 
2.1 Operator product expansion 

In the Bjorken limit, — > cxd, x fixed, the integral in eq.(^ is dominated by the integration region 
near the lightcone 0, because in this region the phase in the exponent in eq.([l|) becomes stationary. 
The external momentum q of the hadronic scattering amplitude is highly virtual, that is to say in 
the unphysical Euclidean region. Thus, we can use dispersion relations together with the operator 
product expansion for a formal expansion of the current product in eq. (||) around the lightcone ~ 
into a series of local composite operators of leading twist. 

The optical theorem relates the hadronic tensor in eq. (I) to the imaginary part of the forward 
scattering amplitude of boson-nucleon scattering, Tf^^, 

Wf,u{p,q) = ^ImTf,u{p,q) ■ (6) 
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The forward Compton amplitude T^u has a time-ordered product of two local currents, to which 
standard perturbation theory applies, 

T^Ap,q) = i/d^^e'''-(P|r(jt(z)J,(0))|P). (7) 

In terms of local operators for a time ordered product of the two electromagnetic or weak hadronic 
currents such as in eq. (0) the OPE reads 

dSe'^-^r(4(z)J,,(0)) = (8) 



N,j 



N 



Qui Qu2\^ ^ r^N 



g2 





{jjuni\9u2pi2l 9uipLiQ.u2Q.ii2 9u2fM2QuiQiJ.i ~^ 9uiu2QiJ.iQii2^ ^2,j 

(7;,3-'?MivO^''^''"-'''^^(/^') + higher twists, 



_,• 1/31/4 ( Q 

'^^UiU2lllUs9 Qu4Qfl2'-'3,j I 2"'*-^ 



where j = Q,q, g. Here, all quantities are assumed to be renormalized, /i being the renormalization 
scale. Higher twist contributions are less singular near the lightcone z-^ ~ and suppressed by powers 
of l/Q^. They are omitted in eq.(|8|). Thus, the sum over N in eq.(|8|) extends to infinity and runs only 
over the standard set of the spin-A^ twist-2 irreducible, symmetrical and traceless flavour non-singlet 
quark operators O", and the singlet quark and gluon operators O*^ and O^. These are defined by, 

Qa,{/.i,-,Miv} = ^aV^iD'^^ ...2)^}^^ a = l,2,...,(nf -1), (9) 

(9q,{m,-,Miv} = ^J^^'^D^'^ ■■■D'^^^i;, (10) 

Qg.f/^l. - .Miv} _ p{uiJ.l JJIJ.2 . . . p/J-Nu} _ ^-l^-j^-^ 

Here, V defines the quark operator and F^'^ the gluon operator. The generators of the flavour group 
SU{nf) are denoted by A°, and the covariant derivative by -D'^. It is understood that the symmetrical 
and traceless part is taken with respect to the indices in curly brackets. 

The spin averaged matrix elements of these operators in eqs.@-(^) sandwiched between some 
hadronic state are given by 

(P|Oi,{w,...,/^iv}|p) = p{Mi...p/^iv}^i^^ 1^^^ ^ (12) 

where hadron mass effects have been neglected. The anomalous dimensions of these operator matrix 
elements in eq.([T^) govern the scale evolution of the structure functions. They are finite quantities as 



well as the coefficient functions C-^j multiplying these operator matrix elements according to eq.(^). 
Both are calculable order by order in perturbative QCD in an expansion in the strong coupling constant 
as- 

The operator matrix elements themselves as given in eq.(|l2|) are not calculable in perturbative 
QCD. However, they can be related to the distributions qi, q^ of quark and anti-quark of fiavour i 
and to the gluon distribution g in the hadron. Let us briefiy note, that the general structure of these 
densities allows for three independently evolving types of non-singlet distributions q^^ and q^g and 
one quark singlet distribution qg. The three non-singlet distributions are the flavour asymmetries 

qns,ij = qi ± qi - (qj ± qj) > (i3) 

and the sum of the valence distributions of all flavours, 

qns = E(q^-q^)' (14) 

i=l 
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while the singlet distribution is simply the sum of the distributions of all flavours, 



i=l 



(15) 



With the operator matrix elements in eq.(12) we can write the OPE of eq.(^) as 
(iSei^-^(P|r(jt(z)J,(0)) |P) = 

m 



(16) 



N 



/X2 



-a 



3,j 



Ai 



+ higher twists, 



which is an expansion in terms of the variable {2p • q)/Q'^ = 1/x for unphysical x — > 00. The final 
connection to DIS structure functions in the physical region < x < 1 is achieved by taking Mellin 
moments of eq.(16) and using the optical theorem eq.(|^), which relates the structure functions to 
invariants Tj^p, i = 2,3, L, of the forward Compton amplitude T^,y as follows, 



Im Ti^p{x,Q'^) 



1 

^imr3,p(x,g2) 



1 

2a; 



Fiix,Q^ 



2,L, 



F3{x,Q'' 



where the Tj^p can be projected in leading twist approximation and D = 4 — 2e dimensions, 

„2 



r2,p(x,Q2) 

T3,p(x,Q2) 



{p-q)' 



■p^p'^T.^ip, 



2e 



2-2e (p-g)2^ ^ 2-2e ^ i M^^v/^^yy, 



1 



(1 - 2e)(2 - 2e) 
Then, with the help of eqs.(|T6|)-(p!^) we find 

1 , 



p-q 



Tf,u{p,q) ■ 



(17) 
(18) 

(19) 
(20) 
(21) 



2,L, 



1 



dxx^-^F3{x,Q'^ 



Ai 



P,N 



(22) 
(23) 



which shows that the Mellin moments of DIS structure functions as defined in eq. (^) can naturally 
be written in the parameters of the OPE eq.(^). 

The derivation of eqs.(p^ and ( p3| ) in the dispersive approach uses symmetry properties of T^jy 
under exchange q — > —q, which is a; — > —x. Therefore, dependent on the process under consideration, 
eqs.(p2[) and ( p3[ ) determine only either the even or the odd Mellin moments of F2, Fi and F^. However, 
all moments in the complex plane are fixed by analytic continuation from either the even or the 
odd Mellin moments. This implies that the x-space result for the physical structure functions in the 
range < a; < 1, can be found by means of an inverse Mellin transformation if the infinite set of either 
even or odd moments is known. 

In the standard case of unpolarized electron-proton scattering in the one-photon exchange ap- 
proximation the even moments of F2 and F^ are fixed in eq. (p2|) , while in the case of electroweak 
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interactions with neutrino-proton scattering eqs.(22) and (p3|) determine the odd moments of F2^~'^^ 
and F^^'^'^^ and the even moments of F^^"*"^^ and F^^~'^ , see ref.[25|. We are going to consider ah 
these cases mentioned in order to extract the complete non-singlet coefficient functions at two loops, 
as will be detailed in section ^. 

The sum in eq.(p^ extends over the flavour non-singlet and singlet quark and gluon contributions. 
Notice however, that in eq.(p3|) the singlet operators O^, do not contribute to F3. This is due to the 



properties of and under a charge conjugation transformation |25|. It can also be understood by 
looking at partonic scattering processes. If p and p are partons and anti-partons, charge conjugation 
implies for the cross-sections 

dp = -dp, (24) 



which gives zero for eigenstates under charge conjugation |15|. 



2.2 Renormalization 

We are interested in the calculation of the scale evolution of the DIS structure functions. To that end, 
a short discussion of the renormalization properties of the operators and the coefficient functions in 
eq.(|^ is in order. 

Let us first recall, that the OPE of eq.(^) is an operator statement and therefore both the coefficient 
functions C/^- and the anomalous dimensions of the operators do not depend on the hadronic states to 
which the OPE is applied. The information on the hadronic target is only contained in the operator 
matrix elements Ap ^ in eq.(|T^). It is therefore standard to consider simpler Green's functions in an 
infrared regulated perturbative expansion with the operators sandwiched between parton states 



(p|OJ-,{Mi,...,/.iv}|p) =^{m...pMiv}^. ^ ^ (25) 




where |p) denotes a spin-averaged parton state, being either a flavour non-singlet or singlet combination 
of quarks and anti-quarks or a gluon. 

As they stand the bare operator matrix elements in eq.(|25|) require renormalization. We choose 



dimensional regularization |2?, ^] in D = 4 — 2e dimensions and define the renormalized operators in 
terms of bare operators as 

0a,hare ^aa Q«,ren (26) 

Qi,bare ^ J- Z^'' O''''^'' , j = q, g, (27) 

k=q_,g 

where Z*^" renormalizes the quark fiavour non-singlet operator 0°. In the fiavour singlet case, eq.(|27|) 
denotes operator mixing under renormalization as and have the same quantum numbers. 
The anomalous dimensions 7 determine the scale dependence of the renormalized operators, 

d 



dln/i' 

and they are defined as 



.^rcn ^ ^O'"'', (28) 



where in the flavour singlet case eqs.(28) and (P9|) are understood as matrix equations, Z representing 
the matrix Z^'^ and 7 the matrix jjk- The general structure of these anomalous dimensions is con- 
strained by charge conjugation invariance and fiavour symmetry. In the case of quarks and anti-quarks 
of flavour i,j they can be split up into a valence and part, 

7q,qi = 7q,q, = ^ijl^q + 7qq , (30) 
7qiq,- = 7q,q, = ^ijlqq + 7qq • (31) 
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As one commonly considers matrix elements of operators corresponding to the quark non-singlet 
and singlet distributions q^sij) Qns qs of eqs.(|l^-(|l^), we remark that their scale evolution is 
governed by the following four linear combinations, 



5.n^^ = 5.,(7^q±7^i) (32) 

^ 7qci^ + ^/7qq^ = 7qq " 7qq + ^/ (7qq " 7qq) , (33) 
qs > 7qq^ + "'/7qq^ = 7qq + 7qq + (7qq + 7qq) • (34) 

In our case of unpolarized lepton-hadron scattering, the even moments of F2 determine the linear 
combinations 7cJ^^ and 7c|^^, while the odd moments of -F3 determine 7^^^ and 7qc^^- The individual 
valence and sea contributions are identified by the flavour structure of the diagrams for the structure 



functions, so that eqs.(32)-(|3^) suffice to determine the anomalous dimensions 7™,7^)7qq and 7^= 



Both, the anomalous dimensions 7 and the renormalization constants Z have a series expansion in 
as- We rewrite eq.(p9|) in D = 4 — 2e dimensions in a form that is suitable for an easy extraction of 
the expansion coefficients, 

7K) = (=11 (ji^Z (o„i)) („„1) , (35) 

where f3{as) denotes the beta function, that determines the renormalization scale dependence of the 
running coupling. The renormalization constants Z obey an expansion in 1/e. Although the Z 
contain poles in e, the anomalous dimensions are always finite as e — 0. Thus, to lowest order in Ug 



the anomalous dimensions 7 in eq.(35) are simply expressed through the residue of Z in 1/e. 

To solve the coupled matrix equation defined by eq.(p9D in the singlet case, one should notice that 
to leading order in the matrix Z^^ is diagonal with zW.qq = = 1 and = Z^o^'S-^ = 0. 

This additional information allows for a unique iterative determination of the anomalous dimensions 
order by order in q^. 

The relation between the bare coupling cxf^'^'^ and the renormalized coupling is given by 

a^-^ = Z^^aT. (36) 
with the renormalization constant Z^^ in the minimal subtraction scheme given by 



e 47r 



(f 








2e 





+ .... (37) 



In D = 4 — 2e dimensions the beta function of the running coupling up to two loops |29| is given by 



11 4 

/3i = ^Cl-ACpTFnf-'^CATFnf. 

Finally, to determine the scale dependence of the coefficient functions it is sufficient to notice, 
that the anomalous dimension of the current in the time-ordered product of the forward Compton 
amplitude T^^ is zero due to current conservation. It implies, that the scale evolution of the coefficient 
functions and the operator matrix elements is governed by the same anomalous dimensions, which 
provides us with the renormalization group equation for the flavour singlet and non-singlet coefficient 
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functions, 



+ /3(a.(/u2)) 



d 



0, 



(39) 



fc=q,g 



5jk - Ijkic^sifJ-'^)) 



0, 



(40) 



In eq.(^) we have adopted the conventional notation to collectively denote the non-singlet anomalous 
dimensions of eqs.(^2|) and (^) with and the coefficient functions with Cj^^^. In particular, the 
Q^-dependence of the coefficient functions does not depend on the index a of the non-singlet operator 
in eq.(P) anymore, see for instance ref. |30|] . 

The actual calculation of the anomalous dimensions as defined in eq.(p9|) and the coefficient func- 
tions C-^j in perturbative QCD proceeds as follows. We introduce partonic invariants Tj^p, i = 2,3, L, 
of the forward partonic Compton amplitude in analogy to eqs.([l9|)-(^Tl). By means of the OPE eq.(^), 
these invariants can be written in terms of renormalized operator matrix elements as 



E 

N 



N 

2x1 ^^'^ 



1 \ iv / f^2 \ / 1 \ In 

J- \ ^N,ns /V \ rynn I \ /ins,ren / P 



^2 T^sj^j ^ns ^ 



<ns,ren i \ , rM 2, 



7i,p(x,Q^,as,e) 



(41) 



(42) 



where we have distinguished the flavour non-singlet from the flavour singlet case, i = 2, 3, L. In the 
former case, we used again the collective notation Z^^ for the various non-singlet renormalization 
constants. The left hand side of eqs.(|4l|) and (42) is renormalized by substituting the bare coupling 
constant in terms of the renormalized one as defined by eq.(B6|), ag = as(/x^/A^). The wave function 



renormalization factors for the external quark and gluon lines are overall factors on both sides of the 
equations and drop out. The terms 0{p^) on the right hand side of eqs.(|4l|) and (^) indicate higher 
twist contributions, which we neglect. 

It is known that the gauge invariant operators and mix under renormalization with unphys- 
ical operators [0, 31]. These are BRST variations of some operators or else vanish by the equations of 
motion. However, matrix elements with physical polarization and on-shell momenta of such unphysical 
operators vanish. That is to say, these unphysical operators do not contribute to quantities related 
to physical 5-matrix elements such as the invariants Tj^p which we are going to calculate. Therefore, 
they are omitted in eq.(^). 

Starting with the partonic invariants Tj^p from eqs.(pT[) and (^2|), the renormalization constants Z 



and the coefficient functions C/^- are calculated using the method of projection developed in ref. [I 



This method consists of applying the following projection operator to both sides of eqs.(41) and 



Vn 



QN 



m 



QpiJ.1 ■ ■ ■ dpi 



PN 



(43) 



p=0 



where q^^^ ■ ■ ■ q>^'^^ is the harmonic, that is to say the symmetrical and traceless part of the tensor 



7^1 . . . a^^N 



On the right hand side of eqs.(pl|) and (42), it is obvious, that the A^-th order differentiation in 
the projection operator Vn singles out precisely the A^-th moment which is the coefficient of l/(2x)^. 
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All other powers of l/(2x) vanish either by differentiation or after nullifying the momentum p. The 
operator Vn does not act on the renormalization constants Z and the coefficient functions on the 
right hand side of eqs.(|4l]) and (^) as they are only functions of N, and e. However, Vn does act 
on the partonic matrix elements A-'^ ^. There, the nullification of p effectively eliminates all diagrams 
containing loops, which become massless tadpole diagrams and are therefore put to zero in dimensional 
regularization. In particular, this removes the operator matrix elements ^g'^" and ^5'^^'^ in eq.(42) 



as they start only at 1-loop level. Hence, in the perturbative expansion of ju- only the tree level 



p,N 

diagrams A^'^^'^'^ survive. Finally, the 0{p'^) terms in eqs.(pl|) and (]4^), which denote higher twist 
contributions, become proportional to the metric tensor after differentiation. They are removed by 
the harmonic tensor q^^^ ■ ■ ■ q^'^^ . 

On the left hand side of eqs.([4l|) and (42), Vn is applied to the integrands of all Feynman diagrams 
contributing to the invariants Tj^p. The momentum p is nullified before taking the limit e ^ 0, so that 
all infrared divergences as p — > are dimensionally regularized for individual diagrams. Effectively 
this reduces the 4-point diagrams that contribute to T^^ to 2-point diagrams with symbolic powers of 
scalar products in the numerator and denominator, which we can solve by means of recursion relations 
as will be detailed in section |^. Notice also, that we apply Vn to the projected partonic invariants 
Ti p rather than to T^jy, as Vn would destroy the tensor structure of T^jj. 



To summarize, we find after application of the projection operator Vn to eqs.(41) and (42) 



rpN, ns ( Q_ 



as,e\ Z^^las,-] ^^'^^'''''(e) 



(44) 
(45) 



where i = 2,3, L and the left hand side is defined as 




VNTi^p{x,Q'^,as,e) 



(46) 



p=0 



Eqs.(|44|) and ( |45| ) are central to our approach of the calculation of the anomalous dimensions 
and coefficient functions via the OPE and dispersion relations. They represent a coupled system of 
equations when both sides are expanded in powers of and e. That is to say, the C/^- are expanded 



in positive powers of e and the Z are expanded in negative powers of e. Explicit solutions to eqs. (|44|) 
and ( p5|) in terms of anomalous dimensions and coefficient functions will be given in section 

At this point a few remarks are in order. First of all, eq.(^) does not provide us with the full 
information about the renormalization constants, because and Z^^ are determined only in the 
order a^. This limitation follows directly from the fact that C/^ only starts from order ag since the 

photon couples directly only to quarks. To extract the 2-loop anomalous dimension 7gq and 7gg^ we 
also calculate Green's functions in which the photon is replaced by an external scalar particle (j) that 



couples directly only to gluons [19|. 

These Green's functions can be expressed in partonic invariants T^.p, for which an OPE similar to 
eq.(|^ exists with the same singlet operators and but with different coefficient functions C^^. 
The important point here is that now starts already at order a^. Hence, the T^^^p, provide us 
with the necessary renormalization constants Z^"^ and Z^ of the singlet operators to two loops. The 
vertices that describe the coupling between the external scalar field (j) and the gluons can be obtained 
by adding the simplest gauge invariant interaction term (pF^^F^'^ to the QCD Lagrangian, where F^^^ 
is the QCD field strength. 
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Repeating the steps that led to eq.(p5|) one finds for these partonic invariants T^^p, 




a 



e. 



(47) 



<p,tree 



p = q,g- 



Beyond tree level we have to take into account the overall renormalization 
(j)po- P-'^ p^^ with the renormalization constant Zp2, 



/ \ bare / \ ren 

(f-^-FI;) = Zp. (f-^'^F^,) + . . . , Zf2 



1 



1 - P{as)/{eas 



of the operator 



(48) 



For the partonic invariants p of the scalar particle cj) this implies an additional overall renormaliza- 
tion 



bare 



(49) 



as indicated on the left hand side in eg. (|47|) . The dots in eq.(48) indicate mixing with unphysical op- 
erators, which give vanishing contributions to on-shell matrix elements with physical spin projections. 
The only physical operator, that mixes with (pF^-^^F^^^ under renormalization is a quark mass term, 
rriqV'V') which vanishes in the limit of massless quarks. 

A second remark concerns the parity- violating structure function F3, which is obtained from the 
time-ordered product of one vector current and one axial vector current A^,. The axial current 
contains a 75 coupling and it is well known that this requires some care in the framework of dimensional 



regularization. We define the axial current by the substitution |3J] 



7m 75 



(50) 



from which F3 is projected according to eqs.(15) and (^Tj) . The sum over dummy Lorentz indices in 
the projection such as in e^^'^'^e^afSx ^^uap defines products of metric tensors, which have to 



be considered as D-dimensional objects [18, 35]. 

The definition eq.(|50|) violates the axial Ward identity which is to be restored by an additional 
renormalization. The necessary renormalization constant Za in the minimal subtraction scheme has 
been given in ref.|ff8[. Up to two loops it is. 



Za 



1 + 



as 
Att 



22 



CaCf - -^Cprif 



(51) 



where the expansion is performed in terms of the renormalized coupling at the scale 

In addition, the treatment of 75 in I? = 4— 2e dimensions introduces an extra finite renormalization 
with Z5 . This is derived in the minimal subtraction scheme from an obvious relation between the vector 
and the axial- vector current. 



(^ms^m) 75 



Z^ Rtj^A 



(52) 



where Rg^^ denotes the i?-operation in the MS-scheme to remove ultraviolet divergencies. It has also 
been given in ref.[|l^ and reads up to two loops. 



Z, = 1-^4C.+ ^ 
47r V 47r 



22Cl-^CACF+'^CFnf 



(53) 
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where again the expansion is performed in the renormahzed couphng as at the scale fi . These 
additional renormalizations in eqs.(51) and ( |53D have to be taken into account when calculating T^^- 
In eq.(H3) one has to substitute on the left hand side for T^'""^ 



(54) 



This concludes our review of the OPE formalism to calculate Mellin moments of DIS structure 
functions. 

3 Methods 



In this section we will discuss the methods used to obtain the Mellin moments of the DIS structure 
functions. 

The idea is to determine reduction identities based on sets of derivative equations for the A'"-th 
Mellin moment of a given diagram in dimensional regularization [p^, 28]. This is done such, that it is 



possible to set up systematically recursion relations in the Mellin moment N. Solving these recursions 
leads to multiple nested sums, which successively can be expressed in terms of harmonic sums, the 
basic functions in Mellin space. This method of recursions was previously used in ref.|jl^ and indeed, 
a number of the relations we give can also be found there. 



3.1 Reduction identities 



Let us begin with a classification of the relevant topologies following the notations of refs.|36 



We will begin with 2-point functions of external momentum q and subsequently dress them up by 
inserting two additional external legs with momentum p. 

The basic topologies of a 2-point function up to two loops are shown below. They are given by the 
basic 1-loop diagram LI and the 2-loop topology named Tl, where the off-shell g-momentum flows 
from right to left through the diagrams. 





Figure 1: The basic basic topologies of a 2-point function up to two loops. Left: The basic 1-loop 
diagram LI. Right: The 2-loop topology Tl. 



The other 2-loop topologies denoted T2 and T3 in refs.[^, 37| are just special cases of Tl, where 
one of the lines 1, . . . , 5 is missing. They reduce to the convolution or to the product of two 1-loop 
integrals. 

For the calculation of the A^-th Mellin moment of deep inelastic structure functions, we need to 
consider 4-point diagrams with external momenta = 0, and q. Then, all topologies up to two loops 
that contribute can be constructed from the diagrams in fig.|l] by attaching two p-dependent external 
legs in all topologically independent ways to the various lines. We call these composite topologies. For 
the calculation we need to construct only a single set of programs that can deal with the A^-th Mellin 
moment of the various composite topologies. 
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Because the method rehes on expressing composite topologies into simpler composite topologies 
we introduce also the notion of basic building blocks. These are composite topologies in which both 
p-lines attach to the same line in the basic topology. At the two loop level there are two basic building 
blocks: 

Till 




1 



1 

{plY{plY{plY{piY{{p+p^Yr ■ 



(55) 
(56) 



In the pictorial representation that we use for the diagrams we indicate the p-flow through the diagrams 
by fat lines. The numbers indicate the number of powers of denominators in the indicated momentum. 
If the same line has two numbers, the second one indicates the number of powers of 2p ■ pi in which i 
is the number of the line in the corresponding basic topology. Hence 




d^'pi d^'pi 



{2p-pi)"' {2p-p5Y 

{p^Y{p^Y{p^Y{plYip^) 



(57) 



We will need the composite topologies Tin, ^"112, TI13, TI14, TI15, Tlig, Tin, Tlss, TI56 and Tlgy. 
All others are related to these by symmetry operations. 

The first observation is that when the external legs 6 and/or 7 are involved, one can expand the 
propagator(s) in p + q with the formula 



{p + qY' 



2n 



E 



2\n 



HYi 



n + i 



l^{2p-qy 



2\i 



(58) 



If we use n powers of 2p ■ q from this expansion, we need only N — n powers oi 2p ■ q from the rest 
of the diagram. Hence such a composite topology becomes a single sum over a simpler topology in 
which p flows through one line fewer. This solves 4 of the 10 composite two loop topologies. 

The way to deal with the other 6 composite topologies is to use all different variations of the 



integration by parts identities and some scaling identities. The integration by parts identities |27 
are of the type 



d''pid''p2 ^ p'; Tl 



d 



kl 



0. 



(59) 



because we deal with a total derivative here. The momenta pi and pj can be equal to any of the 
internal momenta. Additionally pj can be equal to p or q. The scaling identities involve applying one 
of the operators 



P' 



_d_ 

dq^" 



P' 



_d_ 

dpi" 



(60) 



both inside the integral and to the integrated result in Mellin space. The fourth operator of this 
kind cannot be used, because it leads to an inequality when applied to p"^ = 0. Sometimes a third 
type of identities can be useful. It is along the lines of the Passarino-Veltman decomposition in 



formfactors |3£]: 



J d^pid^pa Tiki 



^ -'fcZ ^kl 



(61) 



First the two formfactors I*-'^-' and 1^^^ are determined by contracting eq.(^) either with q^ or p^. 
Then the relevant identity can be obtained by taking the derivative with respect to q^. This can be 
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done because the formfactors are just combinations of powers of 2p ■ q, and a scalar factor. For 
these kind of integrals this method was first introduced in ref.[^]. It is especially useful for the three 
loop integrals. Here we will not need it. 

Using all the above equations, and eliminating all integrals that we do not need, for each topology 
one is left with a solvable equation. Sometimes however some intermediate integrals occur that need 
to be solved by the same methods. First we give the two basic building blocks Tin and Tlssi 



n,N1 n-1,N-1 1 




10 1 1 



V = (n-N-2)/e, 



1 1 





n 




N 


"i 


1 



1 + (-1) 



N 



2(A^+l-2e)(iV-2n+2-2e)' q-q^ 



-(7V-l)(iV-2n+4-4e) 



1 1 



1 1 



-{N-l) 



+8p.p3^(Ar-i) 
q-q 



2p-q 




2 



2 1 



q-q 
q-q 



+-^(4(n-l)(3A^-4n+3-8e)+8iVe-16e 



+Ap-q{2n-N -2+2e) 




{m-An-Ae) 




(63) 



1 1 



Because in both cases the first term on the right hand side is of the same type as the term on the left 
hand side, these equations define proper recursions. All other terms on the right hand sides can be 
computed directly by conventional methods, see for instance refs. 



37| 



Of course there could be powers of 2p-pi present, or different powers of the various denominators, 
but straightforward application of the various relations allows one to reduce all these to the above 
integrals. In the case of Tl^^ one may have to use symmetry considerations. 

The composite topology TI12 is the easiest topology with p-momentum flowing through more than 
one line. Here we start with applying the scaling equation that is based on taking p^^d/dp^ of the 
integral. 



1 1 





(64) 



N is the number of powers of p that needs to be generated in the expansion of the denominators. If 
there are m extra powers of 2p-q it is equal to — m. This last diagram can be simply dealt with. 
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using integration by parts in the left triangle of the diagram. The final result is 



1 1 



(iV + 2) 





3 ^ 1 



1 



+ 



1 



2 1 





1 1 

2 1 




(65) 



and we see three trivial integrals and one integral of the type Tin, albeit with different coefficients 
after the expansion in p. 

For the TI14 topology we have to consider two steps. The first is to reduce the integrals containing 

l/{Pi {p + Pi)'^p\ (p + ^4)^). This is trivial, because we can use the scaling identity based on p^^d/dp^ 

which removes either the 1/p? or the l/pi- After that we can reroute the momentum p and we are 

a 1 



left with a diagram of the type 




in which E is normally 1 and d is 2. We have to be 



careful with the signs if the diagram has to be turned upside down. We are now left with basically 
one diagram with the p momentum flowing through the lines 1 and 5. 

Further reduction can happen in a variety of ways. The important thing to pay attention to is 
that during the reduction, one does not generate simpler integrals with an additional factor of 1/e. In 
that case these simpler integrals would have to be expanded to higher powers in e which may lead to 
relatively difficult sums. Hence we follow a rather elaborate scheme that is free of these problems by 
combining all the equations we can write down for the diagram. It results in: 



1 1 




1 1 



{N + E 



1 1 




1 1 





(66) 



(67) 



(68) 



The second of these equations we will need when we treat the TI13 topology. The third equation leads 
to a recursion which for the fourth term in the right hand side is terminated when enough powers of 
2p-q/q^ have been pulled out. All other terms in the right hand side are of a simpler nature. There is 
no e in any denominator. 
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The derivation of the reduction of the TI15 topology is relatively simple. We use the scaling 
identity based on p^^d/dp^ to obtain 



(iV+2) 



1 ^1 




1 1 




(69) 



The second diagram in the right hand side can be treated with eq.(^). Hence we have to consider 
only the topology — ( )— . The crucial identity for this diagram is obtained in a way similar to 



the derivation of eg. (pq) : i.e. we combine a number of triangle and other relations to form: 



n. m ^ 1 
1 ^1 




1 



N + b + n-m-D 



(n 




n+1 ,m 



n 



1 



n+1. m 

1 

11 11 
n,m 1 n. m 1 n.m- 1 1 n.m- 1 1 




+ 




1 

2 




2 1 





1 1 



1 1 



(70) 



After this there are no more serious problems. The fourth and the sixth diagrams are done by rerouting 
the momentum p through the 4-line and then turning the diagrams upside down. This gives them 
the same topology as the fifth diagram. We can do these diagrams by expanding the propagator and 
applying the tables we constructed for the topology Tin. 

One may note also that because N — m\s always positive or zero, the denominator never becomes 
proportional to e and hence there are no complications with extra singularities in this equation. 

In principle TI13 is the most complicated topology, but because of all the work we have done 
already things turn out to be rather easy. First we notice that the only powers of p in the numerator 
occur in terms of 2p-q. All others can be eliminated leaving terms that could be of a simpler topology. 
Then we write again all the various equations that can be obtained with integration by parts and with 
scaling arguments. In this set of equations we reduce all terms that are not of a simpler topology and 
that contain one of the denominators to a higher power. Eventually we have a rather simple relation 
left: 

1 1 

,2p-q^ 




q-q 




(71) 



In the first term we have an extra power of p and the second term is of a simpler topology. It has 
already been evaluated in eq. (|67|) . 

Actually, this result in eq.(^l|) is a special case of a more general relation. Assume we have a 
general diagram with the structure 




{2p- 



N+n 



CN 



(72) 



in which the blob in the center of the diagram can be anything and the parameter n is determined by 
dimensional considerations. By rerouting the momentum p through the outside line and applying the 
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scaling operator p^d/dp^^ after which the p-flow is restored to how it was in the beginning we obtain 
the relation 




It is clear that in the case of A and B being one, the number of lines that contain the momentum p 
will diminish by one. Hence this relation will be very useful in the future. In the special case of Tlia 
in eq.(|7l|) the first two terms in the right hand side are identical and n = 7—D. 

Many of these equations involve recursion relations and their solution can be written down as a 
single sum over the simpler diagram(s). Each of these sums is a socalled single parameter sum. With 
this we mean that the summand after expansion in e is only a function of the summation parameter 
and the upper limit of the summation, the lower limit being either zero or one. In this hides the power 
of the method, because this avoids multiple sums that cannot be done with current techniques. 



3.2 Solutions 

The equations for the various composite topologies are mostly recursion relations, or in another ter- 
minology: first order difference equations. Suppose we have the equation 

a{N)F{N) -b{N)F{N-l) -G{N) = 0, (74) 

then its solution will be 

F{N) = ^^FiO)+j:%^±l^Gi^). (75) 
Uj=ia{j) i=i Uj=ia{j) 

In the case that the functions a{N) and b{N) can be factorized in linear polynomials in N with the 
coefficient of being one and the other coefficients being integers, the products can be written as 
combinations of F-functions. Because a and b may also depend on the parameter e the F-functions 
should be expanded around e = 0. This will lead to factorials and harmonic sums. The diagram 
indicated by F(0) in eq.(75) can be evaluated with standard MINCER techniques 37]. If the 



function G{N) is expressed as a power series in e with the coefficients being combinations of harmonic 
sums in + m and powers of + m, m being a fixed integer, the sum in eq.(^) can be done and 
F{N) will be a combination of harmonic sums in N + k and powers N -\- k with k being a fixed 
integer. 

Eventually, all solutions to the recursion relations, which were derived above, can be written in 
terms of harmonic sums. Thus, before we continue we will give some attention to these functions that 
describe the solutions in Mellin space. 



3.3 Harmonic sums 

Harmonic sums are the basic functions in Mellin space and much about them can be found in refs.l 



23]. We follow here the notations and the definitions of ref.|22]. Harmonic sums are defined by 



Sm{N) 



N . 

E- 

1=1 * 



N 



S.m{N) = J2 



i-iy 



(76) 
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while higher functions can be defined recursively 

N ^ 

Smi,...,mk{^) = 'J7^Sm2 rrifc ; (77) 

1=1 

S.m,,...,m,iN) = J2^-^Sm„...,md^. (78) 

These functions appear, amoung others, when F-functions are expanded in terms of e. But they also 
show up in sums of the type 

n 1 

EH)Xp;3 = -Si,i,iin), (79) 

1=1 

which are rather common in our calculations. 

The weight of a sum is defined as the sum of the absolute values of all the indices rrij. These sums 
form an algebra. Hence the product of two sums with the same argument can be written as a sum 
of terms, each with a single sum of which the weight is the sum of the weights of the original two 
sums. Many sums involving these harmonic sums can be done by automatic algorithms to any level 
of complexity. For us the important sums are 



E 

1=1 



Sp{i)Sg'{n - i) 



i=l 



^ ' ^ i' 

i=l 



i=l 



in which we interpret an S without indices as 1 and p = pi , . . . , pt, ; q = qi, . . . ,q^. The evaluation 



of these sums is described in |22] and has been programmed in the language of FORM The 
important thing to note is that these sums all evaluate into harmonic sums and denominators of the 
same argument as the harmonic sums. In practise this means that with the sums we run into we will 
remain inside the space of harmonic sums and hence the harmonic sums will span the solution space 
of our integrals in Mellin space. 

3.4 Harmonic Polylogarithms 

In x-space the results will be presented in terms of harmonic polylogarithms |24]. These functions 
are related to the multi-dimensional polylogarithms of ref.[42|. They form a natural basis in the 
space of inverse Mellin transformations of harmonic sums. We will mention here their most important 
properties. For more details the reader should consult ref . p^] . 

The harmonic polylogarithms of weight w and with argument x are identified by a set of mi, ■■■,171^ 
indices which can take one of the three values 0, 1, —1. Harmonic polylogarithms are denoted by 
Hmi,...,mw{x) and explicitly, for lowest weight one defines 

Ho{x) = Inx, (80) 

Hi{x) = / =-ln(l-x), (81) 

H^iix) = r^=ln(l + x). (82) 
Jo 1 + x' 
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For their derivatives, one has 

-^H„,{x) = fra{x) , (83) 

ax 

where again the index m can take the values 0, +1, —1 and the three rational fractions fm{x) are given 

by 

/o(^) = -, /i(x) = -^, /_i(x) = -^. (84) 

X \ — X 1 + X 

In general, the harmonic polylogarithms of weight w are then defined as follows, 

i?mi,...,m„(a:) = ;^ln"'x, if ml,...,m^„ = 0, ...,0, (85) 

while, if mi, ^ 0, 

^^■mi,...,m„(a;) = I dz fmi{z) i^ma ,...,m„ (2) • (86) 

Jo 

To provide a link to the standard literature, we give the results for harmonic polylogarithms up to 
weight three in terms of common polylogarithms in the appendix A. 

Just like the harmonic sums the i?-functions form an algebra. This is to say: products of two H- 
functions, H^^{x)H;^^{x) of weight w and v respectively, where fhuj = mi, ...,m^ and 
can be expressed in a sum of if- function of weight w -\- 

H^Jx)Hn^{x) = ^ J2 ^UM)^ (87) 

in which m^ 1+) fiy represents all mergers of m^ and fiy in which the relative orders of the elements 
of rhyj and n„ are preserved. Hence the sum consists of {v + wy./vlwl terms. This can be shown by 
induction using integration by parts and eq.(86). It should be realized that the rules of this algebra 



are complementary to those of the sums in the sense that the algebra of iZ-functions in x is related 
to the extra algebraic rules for 5-sums in infinity, while the extra algebraic rules for ff-functions in 1 
are related to the algebraic rules for 5-sums in N. 

The harmonic polylogarithms may be expanded in a Taylor series with the expansion coefficients 
being harmonic sums. 



Hrnix) = EEcmn^5n(^), (88) 

-I ^ 
«=1 n 

where a = ±1 and 5^ is a harmonic sum of weight v, ft = ni, ...,n^. In general eq.(^) only holds, if 
m = mi, ...,m^ has no trailing zeroes in the index field. Those correspond to factors ln(x), which do 
not admit a regular Taylor expansion. However, trailing zeroes in the index field can be factored out 
by repeated use of the product identity eg. (|87D , such that eq.(p8[) can safely be applied to the left-over 
i?-function, which by construction does not contain powers of ln(x) anymore. 
Finally, we define the Mellin transformation of regular functions as 

f{N) = V(x) . (89) 







The Mellin transformation of 1/(1 — x) and possible powers of logarithms ln(l — x) are regularized in 
the sense of +-distributions. For this we have to extract first the powers of ln(l — x) which correspond 
to leading indices with the value 1. The extraction is similar to the extraction of trailing zeroes. 
Hence: 



dx x^-i 



Hi^,„^l{x)Hrn{x) 
1 — X 











+ 





x^-^H^{x) - H^il)) ^^3piM , (90) 
/ 1 — x 
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in which m does not have a 1 for its first element. 

Together the above relations allow the construction of the Mellin transform of any i?-function. It 
turns out that there is a one to one relationship between functions of the type Hrfi{x)/{1 it x) with 
m having weight w and S'-sums of weight w + 1. Hence it is possible to construct the inverse Mellin 
transform of any result in space that can be expressed in terms of harmonic sums. The complete 



algorithm has been coded in the language of FORM |41]. 

4 Calculation and Results 

In this section, we will discuss our calculation of the DIS structure functions and list our results 
obtained for the complete set of anomalous dimensions, the flavour non-singlet and singlet quark and 
the gluon coefficient functions up to two loops. 

However, first a few remarks on details are in order. All Feynman diagrams contributing to the 



structure functions have been generated with the help of QGRAF |43] and have been stored in a 



database. For the calculation of the even Mellin moments of the partonic invariants T2^p and 
we used a database identical with the one in ref.[|l^], which contains 425 diagrams up to two loops. 
To obtain the full information about T^^^ and the complete non-singlet contribution to ^2^, we have 
generated a new database of 360 diagrams corresponding to the four different structure functions 
Ff^""^ and Fg-Pi^P. 

Both databases have been checked by calculating some lower fixed Mellin moments in an arbi- 
trary covariant gauge with the MINCER algorithm |^^, keeping the gauge parameter ^ in the gluon 
propagator, that is to say, 

, -9^" + (1 - OQ^q" ^9^) 



— le 

and in the final result all dependence on ^ does cancel. 

The calculation of the partonic invariants T^^, i = 2,3, L and requires to consider only 
external quarks and gluons with physical polarization. In the case of external quarks, the sum over 
all polarizations leads to the projection operator ^ that closes the open string of Dirac matrices 
associated with the external quark line. For external gluons on the other hand, the sum over all 
physical polarizations can be done by contracting the external gluon lines with 

p-q [p-qY 

where p is on-shell, p^ = 0. An alternative approach contracts the external gluon lines only with 
—g°'^. To remove the unphysical contributions, one adds extra classes of diagrams, with external 
ghosts instead of external gluons. The latter approach has been used in ref.[^] while we have checked 
that both methods agree for our calculation of the A^-th Mellin moment up to two loops. 

The explicit calculation of the DIS structure functions has been done with the symbolic ma- 
nipulation program FORM |41|. To that end, all recursion relations given in section ^ have been 



implemented in a program, that reduces the Feynman diagrams of DIS structure functions up to two 
loops to multiple nested sums over the basic building blocks. Subsequently the program calls the 
SUMMER algorithm to solve these nested sums in terms of the basis of harmonic sums. To speed 
up the calculation the results for a large number of basic integrals have been tabulated. 

Let us emphasize that the method of recursion relations as described in section |3| allows for nu- 



merous checks at all stages of the calculation by means of the standard MINCER routine [37|. This is 
from a practical view point by far the most powerful feature of our approach, because the debugging 
of all our programs is extremely simplified. 

A final remark is concerned with the analytic continuation of our A^-space results. In the framework 
of the OPE we calculate either the even or the odd Mellin moments for a given structure function 
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depending on the operators that contribute as detailed in section |2[ This is sufficient to determine 
all moments in the complex plane by analytic continuation. In order to do so we first perform 
the inverse Mellin transformation of our A^-space result to obtain the complete expression in x-space. 
This mapping is unique provided we fix all factors of (—1)'^ according to whether we started from 
even or odd Mellin moments. Subsequently, we may take this x-space result and execute another 
Mellin transformation back to A^-space to obtain the final answer in terms of harmonic sums valid for 
all non-negative integer Mellin moments. This step of restoring all those factors of (—1)^ which are 
due to analyticity is particularly important for the determination of the flavour non-singlet anomalous 
dimensions and coefficient functions at two loops as will be detailed below. 



4.1 Renormalization and mass factorization 

Let us now concentrate on the issue of renormalization and mass factorization for the Mellin moments 
of the DIS structure functions. We wish to show explicitly how to extract the anomalous dimensions 



and coefficient functions from eqs.(44), (|45| ) and (47) if the partonic invariants and are 
expanded in powers of Og and e. In particular, we will obtain the 2-loop anomalous dimensions 7gq 
and 7gg^ from the invariants of the scalar particle. 

The calculation is performed in dimensional regularization D = 4— 2e, [p^, p^. We briefly recall the 



procedure outlined in section y and in refs.[14, 16, 19|. The sum of all Feynman diagrams contributing 



to a given partonic invariant or contains only ultraviolet and collinear divergences. The 
ultraviolet divergences need to be removed by coupling constant renormalization changing the bare 
Q,bare ^j^g renormalized a<j as deflned in eq.(^6|). If necessary, the additional renormalizations 
associated with the axial current, 75 and the interaction term (fyF^^^F^'^ of the scalar particle (j) have 
to be taken into account as described in section § and given in eqs.(p^) and (^4|). 

Then one is left with the collinear divergences associated with the partonic initial states. Those 
need to be removed by mass factorization changing the bare densities of partons in the hadron to 
renormalized ones. This is the same as renormalizing the operator matrix elements Ap ^ as done in 



eqs.(y), (|5D and 

In the following, we list the results for the Mellin moments of the parton invariants order by order 
in as- To that end it is useful to deflne the following expansions in powers of as for the coefficient 
functions and the anomalous dimensions, 



n=0 



n=0 



47r / 

asV" (n) 



C^A-s) = E C\ (94) 



and, since the operator matrix elements ^p'^^*^ factorize after application of the projector Vn, also 
for the parton invariants, 

= {t^!S+tS+T^S + ...)aI%-, (95) 

where all left-over singularities in the T-'-p are of collinear nature. The procedure of mass factorization 
works iteratively order by order, both in as and in the regularization parameter e, such that all 
physical quantities can be uniquely extracted. Notice in particular, that the gluon tree level operator 



matrix element ^g'^'^*' is (1 — e)const7v. If factorized as in eq.(95), the factor (1 — e) due to the gluon 
polarization in D dimensions is essential for a proper determination of the 1-loop coefficient function 
C2^g. As discussed in ref.|16], the omission of the factor (1 — e) accounts for some of the differences 



with ref.jT^ for the structure function -Fl.Q 
^ The discrepancies were corrected in ref . . 
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All results are given in dimensional regularization in the MS-scheme with as being the renormalized 
quantity according to eq.(pl). At leading order, we have normalized, 



.(0) 



(0) 



1 



.(0) 



.(0) 



(0) 



2,q 3,q ' 2,g L,q L,g 

At first order in as, we have to expand up to order e and find. 



0. 



(96) 



(1) 



2,q 



(1) 



2,g 

r(i) 

-L,q 

r(l) 



as c- 
An 



as n 
An 



47r 



g 7qq ''2,q ^ *="2,q 



g 7qq ^3,q ^ *="3,q 



^7qg +C2^g + ea2^g 



(97) 
(98) 
(99) 
(100) 
(101) 



where the factor is defined by 



= exp (e{ln(47r) - 7e}) 



(102) 



At second order in Og, we need to split up the contributions into flavour non-singlet and singlet 
parts. Allowing for electroweak interactions, one can consider in the non-singlet case the structure 



functions of four different physical processes, F2 
distinguish even and odd moments. We have 



uP±uP 



and R 



uP±uP 



This implies that we also need to 



p(2),ns,± 
2,q 



T, 



(2),ns,± 



3,q 



p(2),ns 
-L,q 



+- 



1 n 



27qq 



(i),±,v 



+ 7, 



(0)J1) 
qq 2,q 



(2),ns,+ (2),ns,- (0) (1) 
+ '^2,q =■= ^2.0 -t- 7qq "r 



as 
Att 

1 n 



S 



2e r- 



1 n 



^(0) 
/qq 



-2,q 



/qq "'2,q 



+-^^7S'^'^ + 7(°JcWk 



47r 



7qq '-3,q 
2e 



(2),ns,+ , (2),ns, 



-3,q 



3,q 



^ 7qq "3,q 



\ /qq ^L,q / 



+ C 



(2).ns^ .0)^(1) 



L,q 



qq "L.q 



(103) 



(104) 



(105) 



where the it-sign denotes even or odd moments in the expressions for T2 ^''"'^'^ and T^^^'^ . We want 
to emphasize this distinction for even and odd moments as done in eqs.( |l03|) and ( |104| ) is due to the 
fact that we deal with different scattering processes. 

It is relevant for the 2-loop anomalous dimensions entering in eqs.( |103| ) and ( |104D because of flavour 
symmetry breaking in the anti-quark distributions. We use the notation of eq.(|3^) 



/qq 



'qq 



±7. 



(i),v 
qq 



(106) 
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It is also relevant for the 2-loop coefficient functions, which decompose for even and odd moments as 



.(2),ns 
„(2),ns 



„(2),nsH 

„(2),nsH 
-2,q 



+ -irc: 



AT (2),ns,- 



2,q 



AT (2),ns,- 



3,q 



(107) 

(108) 



After analytical continuation of our A^-space results we can reconstruct these physical signs (—1)^ in 
eqs.(106)-(108). That is to say, we can then determine 7^ and 7^ by taking the sum or the difference 
of 7^'^ and 7qq'^- Analogously, the individual contributions to the non-singlet coefficient functions at 
two loops, Cg^q"*", C2q~, c^^^ j a^d Cg^q" , are obtained, but again the reconstruction works only after 
analytical continuation. 

Note also, that in the equations above the partonic invariants Tg ^ and Tg q^'"'' are always understood 
to be renormalized according to eq.(54). 

(2) s 

We decompose the quark singlet invariant T- q into non-singlet and pure-singlet contributions. 



p(2),s 
-2,q 



m(2),nS,+ rTn(2),pS 

-'2,q "'"-'2,q ' 



rTn(2),S _ m(2),nS m(2),pS 



(109) 



and for the pure-singlet contributions we have at second order in as, 



r(2),ps 

-2,q 



71/ 



4-K J 



2e . 



1 n 



e2 12 



(0)^(0) 



I'^qg '''gq 



+ 



1 /1^(1)>+,S , ^(0) (1) I , (2),ps ^(0) (1) 
^ i r,7qq + 7gq C2,g f + (^2,ci + 7gq «2,g 



e 12 



r(2),ps 



Uf 



2e 



\ 'gq ^L,g / 



+ C 



(2),ps (0) (1) 
I Ig-Q ^ 



X,q 



gq L,g 



(110) 



(111) 



In eq. (|110| ) we used again the notation for the anomalous dimensions of eqs.(p3[) and 



Iqq 



'qq 'qq 



(112) 



The pure singlet contribution to T3 q^''^'' vanishes, implying that we have 



^(i),-,s 
'qq 



(1),S (i),s 

7qq -7qq'' 



0. 



(113) 



For the gluonic invariants, we find at second order in Os, 

2e 



n{2) 
-2,g 



^^^2 /^if^y [1 fl (0) f (0) ^ {0)\ _ 1^ (0) 

4vr/ ' J |.e2 I2 '^''s V '^i'^ ^ J 2' ^"^^ 



1/1^(1) +^(0)J1) , ^(0)J1)1 , J2) , (0) (1) , (0) (1)' 
^ I 2'qg ^ 7gg '-2,g ^ 7qg ''2,q ( ^ ^2.e: ^ I Re. "2.2 ^ Iqk ":~ 



r(2) 



2e 



47r / 



2,g /gg "'2,g -r rqg "2,q 

1/^(0) Jl) +^(0)c«\ , J2) , ^(0) (1) , ^(0)J1)' 
f t^gg ^L,g + 7qg Ci,q I + C^,g + 7gg O^^g + 7qg «L,q 



(114) 



• (115) 



Finally, we can give the expressions for the partonic invariants of the scalar particle (p. We find at 
leading order. 



J.(0) _ Q 



.(0) 



1 



(116) 
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In the equations below the partonic invariants T^^p and T^^p are always understood to be renormalized 



according to eq. ([49[) . Then we obtain at first order in a^, 



.(1) 



.(1) 



4tt [ Q2 



^(0) + ^(1) 



_ /gg -r t-j, „ 



where we only have to expand up the finte terms in e. 
Finally, at second order in a^, 



.(2) 



2e . 



An) [q\ 

1 n 

e 



1/1^(0) Uo) +^(0) 



^ ^ 2 ^'^ ^.g ^ 'qq (/-.q 



.(2) 



+ 



47r 

1 n 



2e . 



1 n 



r^(o)^% 1 (0) (o)_i^ (0) 



(117) 
(118) 



(119) 



(120) 



This concludes our brief discussions on the extraction of the anomalous dimensions and coefficient 
functions from eqs.(||), (|5|) and (^ by means of eqs.(|96|)-( |120D . 



4.2 Results in Mellin space 

We are now ready to present our results, which provide all necessary ingredients for the solution of 
the renormalization group equations ( |39|) and (^) for the flavour singlet and non-singlet coefficient 
functions as detailed for example in ref.|^5|. However, as is well known, a complete solution to NNLO 
for the scale evolution of the coefficient functions requires also the still unknown anomalous dimensions 
7pp'' at 3- loops. 

To summarize, our results for the anomalous dimensions agree with the ones published in refs.|Q, ^, 
@) |5|, ^, 1^, ^] while our results for the coefficient functions agree with those of refs.[Q, |l^, [ll, [T5|, [l^ . As 
far as the earlier calculations of the longitudinal structure function at two loops of refs7p|JlO|, |11|, p^ j 
are concerned, let us mention that there is only complete agreement for c^^ and with ref.[^]. 
Part of the discrepancies with ref.[^] have been explained above. 

For the anomalous dimensions, we find at 1-loop, 



7S(iV) 
7£)(iV) 



Cf{-3 + 2Si{N-1) + 2Si{N+1)} , (121) 

nf{2Si{N-l)+8Si{N+l)-ASi{N+2)-6Si{N)} , (122) 

Cf {45i(iV-2) - 8Si{N-l) - 2Si(iV+l) + 6Si(iV)} , (123) 

CA{^SiiN-2) - 8Si{N-l) - 8Si{N+l) + 4Si{N+2) + 12SiiN)} - Po . (124) 



At two loops, we obtain. 



7i?^(iV) = (-1)^ X 



CpCa|(-1)^ (- 



(125) 



17 /268 

h 4C3 + 

6 V 9 



8C2)Si{N-l) - ^S2{N-1) + 8S3{N-1 



-4S-3{N-l) - 4S^3{N+1) + 8S-3{N) + ^S-2{N-1) + y 5_2(A^+1) 
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-fs.2{N) + (1^ + 4C2)s_i(iV-l) - -4C2)5_i(iV+l) 
+ (T-8C.)^-i(A^)) 



4 44 40 

--S_i(iV-l) + -g-S_i(iV+l) - -jS^iiN) 



+ C|.|(-l)^ (-^ - 8C3 + 165i(A^-l)C2 - 16Si,2(A^-l) + 12S2{N-1) - 16S2,i{N-l)) 

-4S^3iN-l) - 4S-3iN+l) + 8S^3iN) + 85_2(iV+l) - 85_2(iV) + 8S.2,iiN-l) 
+85_2,i(A^+l) - 165_2,i(iV) - (20 + 8C2)5_i(iV-l) + (20 - 8(2)5.1 (A^+1) 

+165_i(7V)C2 + 8S_i,2(iV-l) + 8S_i,2(A^+l) - 16S_i,2(iV) 

Tg''' W = (-1)^ X (126) 
Cf(Cf - ^) jsCs + 85_3(A^-1) + 8S-3{N+1) - 8S-2{N-1) - 85_2(iV+l) 
+165_2(A^) + 165_i(iV-l) - 165_i(A^+l) - 85i(iV-l)C2 - 85i(iV+l)C2 
-165i,_2(A^-l) - 165i,_2(A^+l)| , 

7ff''(iV) = 7S''W = (-l)'^x (127) 

f 46 
C^n/^ -45_3(iV-l) - 45_3(iV+l) + 8S_3(A^) - 25_2(A^-1) - y5_2(A^+l) 

+y5_2(A^+2) + 125_2(iV) - ^5_i(iV-2) + ^5_i(iV-l) - ^5_i(iV+l) 
+^S_i(Ar+2)-8S_i(Ar)|, 

7£)(iV) = (-l)^x (128) 

CFnfi^4S-s{N-l) +8S-3{N+l) - 165_3(iV+2) + 45_3(7V) - 65_2(Ar-l) 

-85_2(iV+l) + 165_2(A^+2) - 2S-2{N) - 8S-2,iiN~l) + 165_2,i(A^+2) 
-85-2,1 (A^) + 285-1 (Af-1) - 185_i(Af+l) -405_i(Af+2) +305_i(Af) 
-165_i,i(Ar+2) + 165_i,i(iV) + 85_i,i,i(iV-l) - 165_i,i,i(A^+2) + 85_i,i,i(A^) 

-85_i,2(7V-l) + 165_i,2(iV+2) - 85_i,2(7V)| 

f 272 
+ CAUfl -85_3(7V-1) - 165_3(7V+1) + 245_3(A^) - 45_2(7V-1) - — 5_2(iV+l) 

+^5_2(iV+2) + 365_2(iV) - ^5_i(iV-2) + + 4(2)5.1 (iV-1) 
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+y5_i(Af+l) + - 8C2)5-i(iV+2) - (92 - 4C2)5_i(iV) + 165_i,i(iV+2) 
-165_i,i(iV) - 85_i,i,i(iV-l) + 165_i,i,i(iV+2) - 85_i,i,i(Af) - 45i(iV-l)C2 
-165i(Ar+l)C2 + 85i(A^+2)C2 + 125i(iV)C2 - 85i,_2(A^-l) - 32Si,_2(A^+l) 

+165i _2(A^+2) + 245i _2(A^) 
7«(iV) = (-l)^x (129) 

{92 
165_3(iV-l) + 85_3(iV+l) - 245_3(iV) + 485_2(iV-l) + -jS-2{N+l) 

32 

-— 5_2(iV+2) - 68S_2(Ar) + 165_2,i(iV-2) - 85_2,i(Ar+l) - 85_2,i(iV) 

-(4 + 8C2)5_i(A^-2) + li^5_i(Ar-l) + (36 + 4C2)5_i(iV+l) - lI^5_i(iV+2) 

-(^ -4C2)s_i(iV) + ^5_i,i(iV-2) - ^5_i,i(Ar+l) - y5_i,i(iV) 
-165_i,i,i(7V-2) + 85_i,i,i(Ar+l) + 8S_i,i,i(A^) + 165_i,2(A^-2) - 85_i,2(iV+l) 
-85_i,2(A^) - 85i(iV-2)C2 + 165i(7V-l)C2 + 45i(7V+l)C2 - 125i(iV)C2 

- 165i,_2 (iV - 2) + 325i,_2 (iV - 1) + 85i,_2 (iV + 1) - 245i,_2 (iV) I 
+ CFnfl^^S.i{N-2) - ^5_i(iV+l) - ^5_i(iV) - ^S.i,,{N-2) + ^5_i,i(iV+l) 
+^5_i,i(iV)| 

+ C|,|-85_3(iV-l) +45_3(7V+1) + 45_3(7V) -8S_2(iV-l) - 14S_2(A^+1) + 225_2(iV) 
-105_i(Af-l) - 145_i(iV+l) + 245_i(iV) - 245_i,i(iV-2) + 205_i,i(A^+l) 
+4S_i,i(7V) + 16S_i,i,i(iV-2) - 8S_i,i,i(iV+l) - 85_i,i,i(iV) 

7«(iV) = (-l)^x (130) 
CFnfl^2{-l)^ - 85_3(iV-l) - 85_3(iV+l) + 165_3(^) + 125-2(A^-1) + 205_2(A^+1) 

-325_2(iV) - ^-S_i(iV-2) - ^S_i(iV-l) + ^S_i(iV+l) - ^5_i(7V+2) 

+165_i(A^)| 

+ CAnfl^{-lf(l - ^S,{N-l)) + ^5_2(iV-l) + ^5_2(iV+l) - y5_2(iV) 

92 R Ifi 99 4n 1 

+-5_i(iV-2) + -S.,{N-1) + -S_i(iV+l) - y5_i(iV+2) - y^-iWj 

+ C'l|(-l)'^(-y - 4C3 + + 8(2)^1 (iV-1) - 165i,2(iV-l) - 1652,i(iV-l) 

+853(7V-1)) + 4C3 + 8S_3(iV-l) + 16S_3(iV+l) + 165_3(iV+2) - 32S_3(7V) 

24 



+^5_2(iV-l) + 445-2(iV+l) - ^5-2(iV+2) - ^5-2(iV) + 165_2,i(iV-2) 
+165_2,i(iV-l) - 165_2,i(A^+2) - 165_2,i(A^) - 85_i(iV-2)C2 
-(^ + 8C2)5-i(iV-l) - ^5-i(iV+l) + 8S_i(iV+2)C2 + + 8C2)S-i{N) 
+165_i,2(iV-2) + 165_i,2(A^-l) - 165-i,2(iV+2) - 165„i,2(iV) - 85i(iV-2)C2 
+165i(A^-l)C2 + 165i(Af+l)C2 - 85i(iV+2)C2 - 245i(iV)C2 - 165i,_2(A^-2) 

+32Si,_2(iV-l) + 325i,_2(iV+l) - 165i,_2(iV+2) - 485i,_2(iV) 



For the coefficient functions, we obtain at tree level, 



At 1-loop we find, 



.(0) 
-3,q 



(0) 



(131) 



(132) 



.(1) 
.(1) 



(iV) 
(iV) 



(AT) 
(A) 



(133) 



{-9 - 35i(A^-l) + 45i(A+l) + 25i(A) + 25i,i(A-l) 
+25i,i(A^+l) - 252(A^-1) - 252(A^+1)} , 
c^lliN) + CF{2Si{N-l)-2Si{N+l)} , 

nj {2^1 (A-1) + 325i(A+l) - 165i(A+2) - 185i(A^) + 25i,i(A-l) (134) 
+85i,i(A+l) - 45i,i(A^+2) - 65i,i(A) - 252(A^-1) - 852(A^+1) 
+452(A+2) + 652(A)} , 
4Ci.{5i(A^+l) -5i(A^)} , (135) 
8nf {25i(A^+l) - 5i(A^+2) - Si{N)} . (136) 

(2),ns (2),ns (2),ps (2) {2),ns 



Our results for the Mellin moments of 2-loop coefficient functions C2 q ,03 q , C2 q , C2_g q 



„(2),ps 
-L,q 



and c^^g can be found in appendix B. 



4.3 Results in x-space 

Let us now present our results in momentum space. The scale evolution of the DIS structure functions 



in x-space has been discussed in ref.|46]. It is governed by the Altarelli-Parisi splitting functions, 
defined as 



(137) 



where the conventional relation between the anomalous dimensions and the splitting functions in 



eg. ( 137 ) involves a relative sign. In the following, all expressions which diverge for x — > 1 are under- 



stood in the sense of +-distributions as defined in eq.(90). 
For the splitting functions, we find at 1-loop, 



CF{2pqq(x)+35(l-x)} 

2n/pqg(x) , 
2CFPgg(,x) , 

4CAPggix) + Pq5{1 - x) , 



(138) 
(139) 
(140) 

(141) 
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where we have introduced the following polynomials, 



Pqq(x) = (142) 
1 — X 



) = 1-23; + 23;% (143) 
2 

Pgq(x) = - - 2 + X, (144) 

X 

Pgg{x) = -^— + --2 + x-x'^. (145) 
^ 1 - X X 

At two loops, we obtain, ^ 
P«'V(x) = (146) 

- ^) - ^(^ + ^0(x))pqq(x) - + ^(2)S{1 - x)^ 

+ CpC^jy (1 -x)+ 4(1 + x)/7o(a:) + 4(^ - C2 + y^o(x) + /7o,o(x))pqq(x) 

+ (f + |C2-12C3)5(1-X)} 
+ C||-20(l - x) - 2(3 + 7x)Ho{x) - 8[^Ho{x) - Fi,o(x) - H2{x))p^^{x) 

-4(1 + x)Ho,o{x) + - I2C2 + 24(3)5(1 - x)| , 

Pii'^ix) = (147) 
Ci.(Cir - ^) 1 16(1 -x)+ 8(1 + x)Ho{x) + 8(/7o,o(x) - 2F_i,o(x) - C2)pqq(-x)| , 

= Pii'^'ix) = (148) 
njC^|-4 + + 12x - ^x^ - 4(1 + x)Fo,o(x) + 2(1 + 5x + ^x2)/7o(x)| , 

P(i)(x) = (149) 
n/CF|8 - 18x - 2(1 - 4x)Fo(2;) - 8i?i(x) - 4(1 - 2x)Hofi{x) 

+8(^ - C2 + Ho{x) + Ho,o{x) + Hi{x) + /i'i,o(2;) + Hi^i{x) + H2(x))pqg(x)| 
^ [364 80 28 „/109 11 , , , , , , ,\ , , 



(19 - 68x)/7o(^) + 8Hi{x) - 8(1 + 2x)Fo,o(x) - 8//_i,o(x)pqg(-x) - I6XC2 



^The definition of 5*2 in ref.|6| contains a typographical error. The lower integration boundary {1 + x)/x should read 
x/{l + x), see for instance eq.(61) in ref.|^. Also the timelike quark-quark splitting function in eq.(12) of ref.^ contains 
a typographical mistake. The term (10 — 18x — 16/3x^) ln{x) should read (—10 — 18x — 16/3x^) ln(a:). 
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pW(x) = (150) 
nfCFl^~x-^(^-Hi{x))pgq{x)^ 

{112 130 176 1 11 

— + —X + —x^ + 8(- - —Hi{x) + i?i,o(x) + + H2{x))pgq{x) 

-4(12 + 5x + ^x'^)Ho{x) - 8xHi{x) + 8(2 + x)Ho,o{x) - 8iI_i,o(x)pgq(-x) + 16C2| 

+ C|^|-10 - 14a; - 4(2 - x)Ho^o{x) + 2(4 + 7x)Ho{x) + 8xHi{x) 

+8(^i/i(x)-i/i,i(x))pgq(x)|, 

Pg\x) = (151) 
UfCpl -32 + ^ + 16x + ^x'^ - 8(1 + x)Hofi{x) - 4(3 + 5x)Hq{x) - 2(5(1 - x) I 

I oX o I 

+ nfCAl^4{l -x)-f{^-x')- ^pgg(x) - l{l+x)H,{x) - ^<5(1 - x)| 

{268 1 4 
54(1 - x) - — (- - - -(25 - llx + A4x^)Ho{x) + 32(1 + x)Ho,o{x) 

+8(^ - C2 + ^fo,o(a;) + 2i7i,o(x) + 2i/2(a;))pgg(x) 

+8(/fo,o(^) - 2iJ-i,o(x) - C2)pgg(-x) + (f + 12(3)^(1 - x)| . 



1 

;,p(iV) = J dxx^-'ci,p{x). (152) 



The x-space expressions for the coefficient functions are commonly defined by, 

1 



We obtain at tree level, 

= cflix) = S{1 - x) , cflix) = cflix) = cflix) = , (153) 

At 1-Ioop we find, 

4;)(x) = Cp|^ + |x-2(^ + ifo(x) + i?i(x))pqq(x)-(9 + 4C2)5(l-x)|, (154) 

c'ilix) = c^2,1{x)-2Cf{1 + x), (155) 

c2(x) = n/je- 2(4 + i7o(x)+i7i(x))pqg(x)|, (156) 

cS(^) = ^Cfx, (157) 
cW(x) = 8n/x(l-x). (158) 

5sults for the 2-Ioop 
been deferred to appendix C. 



The results for the 2- loop coefficient functions C2^q'^'^,C3^q°^, Cg^q^^i q"''' ^^q''** c^g have 
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5 Conclusions 



In the present paper we have calculated the Mellin moments of the perturbative QCD corrections up 
to second order for the DIS structure functions F2 , -F3 and Fl in leading twist approximation using the 
method of projection [^]. We have presented the analytic results for the 1- and 2-loop contributions 
to the anomalous dimensions of the singlet and non-singlet operator matrix elements and the 1- and 
2-loop coefficient functions of F^ and F^. Our results are in agreement with the literature as has 
been detailed in section ^. 

Our choice to work in Mellin space enabled us to solve the 2-loop integrals in dimensional regu- 
larization, D = 4 — 2e, by means of recursion relations in the Mellin moment N . This approach has 
systematically extended previous work in this direction jl^ and relied on the improved understanding 
of harmonic sums @, Progress was possible in particular due to the development of new 

algorithms for a large class of series that involve harmonic sums. 

The method as we applied it turned out to be very flexible and in the expansion in e it is in 
principle not limited to a certain order. This allowed us to calculate anomalous dimensions and 
coefficient functions at the same time. As a byproduct, we have also performed the calculation of 
all 0(e^)-terms at 1-loops and all 0(e)-terms at two loops, for all unpolarized structure functions, as 
needed to extract the coefficient functions at 3-loops after mass factorization. However these results 
will be published elsewhere [47|. 

We have been able to perform the inverse Mellin transformation of our A^-space results analyt- 
ically and we have given the corresponding expressions in x-space as well. Our x-space results are 
expressed in terms of harmonic polylogarithms, which we believe to be the natural class of functions 
for calculations of DIS structure functions. Again, this has been achieved due to deeper insight gained 
into the subtle interplay between harmonic sums and harmonic polylogarithms [24|. 

We are very confident, that the program presented in this work can be successfully applied to the 
ultimate goal, the calculation of the anomalous dimensions and the coefficient functions of the DIS 
structure functions at 3-loops. Finally we would like to note that a calculation of the perturbative QCD 
corrections to the polarized structure function gi up to second order and in leading twist approximation 
is in progress |47|. 
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Appendix A 

Here we give the results for the harmonic polylogarithms of weight two and three in terms of standard 
polylogarithms , 

Li2(x) = - / — ln(l - z) , Lhix) = / —U2{z) , (159) 

Jo z Jo z 

where Li2(x) is Euler's dilogarithm and Li3(x) the usual trilogarithm, see ref.|4^. For It is obvious, 
that we only need to consider a limited subset at weight two and three. We find at weight two. 



z 

H-ifl{x) = lnxln(l -|- x) -|- Li2(— ; 



(160) 
(161) 
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= -k2 + ^ln^2 
z z 


-ln(l+x)ln2 + Li2f^) , 


(162) 


Ho,-i{x) 


= -U2{-x); , 




(163) 


Hofi{x) 


= - In^ X , 
2 




(164) 


Ho,i{x) 


= l^hix) , 


/1 + t\ 

ln(l - z) ln(l + x) + ln(l + x) In 2 - Li2 , 

V 2 y 


(165) 


Hi,-i{x) 


z z 


(166) 


Hifl{x) 


= — lnxln(l — x) 


- Li2(x) , 


(167) 


Hi,i{x) 






(168) 



At weight three we obtain 
1 



1,-1,-1 



(x) 



i?_i_i,o(x) 
H-,,-,,i{x) 



H-ifl,o{x) 

H-ifl,i{x) 



H-i,i,o{x) 



H-i,i,i{x) 



^\v?{l+x), 
Ca - ln(l+x)C2 + \ ln3(l+x) - Lig 

D 



1+X 



ln2 - ^Cs - ^ ln=^ 2 - i ln(l+x)C2 + \ ln(l+x) In^ 2 



1 



ln2(l+x) ln2 + Li3l 



1+x 



ii'_i,0 -i(x) = -2C3 + 21n(l+x)C2 - 3 ln3(l+x) + ln(l+x)Li2(-x) + lnxln2(l+x) 



+2Li3 



1+x 



- In^ xln(l+x) + lnxLi2(— x) — Li3(— x) , 

^Cs - \ ln(l-x)C2 - ln(l-x)Li2(-x) + ]- ln(l+x)C2 + \ ln3(l + x) 

4 Z ZD 



-lnxln(l-x)ln(l+x) -Li3(l-x) - Li3^--^-j-^^ - Li3^ 



1 



1+x 



+Li3 



l-x\ 



iJ_i,i,_i(x) = -C2ln2+^C3 + ^ln^2 + iln(l+x)C2-^ln(l+x)ln22 



+ ln(l+x)Li2 



1+x 



-2Li3 



l+x\ 



(169) 
(170) 
(171) 

(172) 

(173) 

(174) 



(175) 



-iC2 ln2 - + ^ In^ 2 + ^ ln(l+x)C2 - \ ln(l+x) In^ 2 
+ - hi2(l+x) ln2 - -ln^(l+x) - -lnxC2 + - lnxln^2 - lnxln(l+x) ln2 



(176) 



+^lnxln^(l+x) +lnxLi2^ 



1+x 



Li: 



. (l+x\ 



+ L„(-x) + Li3(^) 



+Li3 



2x 
1 + x 



- Li3(x) , 



iC2 ln2 - - ^ In^ 2 + ln(l-x)C2 - \ In(l-x) In^ 2 

+ ^ ln2(l-x) ln2 - ^ ln2(l-x) ln(l+x) + ^ In(l-x) ln2(l+x) 

- ln(l-x)Li2(^^^ - ln(l+x)C2 + ^ ln(l+x) In^ 2 - ^ ln\l+i 



(177) 
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+Li3(^)+Li3^-i^) 



Ho,-i,-iix) = C3-ln(l+x)C2 + 7ln2(l+a;)-ln(l+a;)Li2(-x)-^lnxln2(l+x) (178) 

D 2 



-Lial 



1 



1+x 



Ho-i,o{x) = -lna;Li2(-.T) + 2Li3(-x), (179) 
Ho,-i,i{x) = ^C2ln2-^C3-^ln'2+iln(l-a;)C2 + ln(l-x)Li2(-x) (180) 

-2 ln(l+a;)C2 + ^ ln(l+x) In^ 2 - ^ ln^(l+x) In 2 + ^ ln^{l+x) 
2 2 6 

+ lnxln(l-a;)ln(l+a;) - - lnxln^(l+x) + W j +Li3(l-x) 

-Li3(-.) + Li3(-^) - L.3(^) - L.3(^) + Li3(x) , 
Ho,o,-i{x) = -Lhi-x), (181) 

HoM^) = ^In^a;, (182) 
Ho,o,i{x) = Usix), (183) 
Hq,i,-i{x) = -lc2ln2-^C3 + ^ln^2+^ln(l+x)C2-^ln(l+ar)ln2 2 (184) 

+ -In2(l+x)ln2 - -ln^(l+x) + ln(l+x)Li2(x) + - lnxln^(l+a;) 

-Us[^) + Li3(-X) + Li3(^) + Li3(^) - Li3(x) , 
Ho,i,o{x) = lna;Li2(a;) -2Li3(a;), (185) 

Ho,i,i{x) = C3 + ln(l-x)C2-ln(l-x)Li2(x)-^lnxln2(l-x)-Li3(l-x), (186) 

iJi,_i,_i(x) = ^C2ln2-^C3-^ln'2-iln(l-a;)ln2(l+a;) + ^ln2(l+x)ln2 (187) 



Hi,-i,o{x) = lc2ln2-^C3-^ln32 + ^ln(l-a;)C2-21n(l+a;)C2 (188) 

+ -ln(l + x) 111^2 - -ln2(l+x) ln2 + - ln2(l+x) + -lnxC2 - -lnxln^2 
2 2 6 2 2 

1 / l+x\ / 1+x 

+ lnxln(l+a;)ln2 - -lnxln2(l+x) -InxW j +Li3( 

+Li3(l-x) - Li3(-x) + Li3(-^) - Li3(i^) - Li3(^) + L.3(X) , 

Hi,-i,i{x) = -C2ln2 + ic3 + ^ln=^2-^ln(l-x)C2 + ^ln(l-x)ln22 (189) 
- In2(l-x) ln2 + ln2(l-x) ln(l+x) + In(l-x) ln(l+x) ln2 
-ln(l-x)ln2(l+x) +ln(l-x)Li2(^^^^ + 21n(l+x)C2 - ln(l+x) In^ 2 

+iln3a+.)-2L,(i±i:)-2L,( iz£), 

i?i,o,-i(x) = ic3-^ln(l-x)C2 + ^ln(l+x)C2 + ^ln-^(l+x)-ln(l+x)Li2(x) (190) 



30 



1 — X 

lnxln(l — x) ln(l+a;) — Li3(l — x) — Lis! ) — Lis 



1+x J \l+x 



+Li3 
1 



1 — x 
1+x 

2 



Hifi^o{x) = --ln^xln(l-x) - lnxLi2(x) +Li3(x) (191) 
= Cs - S„,p(l - x) , 

Hi,o,i{x) = -2C3-21n(l-2;)C2 + ln(l-x)Li2(x)+lnxln2(l-2;)+2Li3(l-x), (192) 

Hi,i,-i{x) = ic2ln2-^C3-^ln'2+iln(l-x)C2 + Jln2(l-x)ln2 (193) 
2 o 2 2 

- In(l-x) ln(l + x) In 2 + ^ In(l-x) ln^(l+x) - ln(l+x)C2 + ^ ln(l+x) In^ 2 
-ll„3(l+.) + Li3(l±H)+Li3(-i^). 

Hi,i,o{x) = C3 + ln(l-x)C2-Li3(l-x), (194) 
Hi,i,i{x) = -^ln3(l-x). (195) 

D 

The function S„_p in eqs.( |178[ ) and ( |191] ) denote the Nielsen functions |4S], defined as 

1 1 

(_l\p+n~l I- J 

S„,p(x) = \/ / - ln"-i(z) ln^(l - xz) . (196) 



pi (n — 1)1 J z 


Appendix B 

Here we present the formulae for the Mellin moments of the 2-loop coefficient functions. We obtain, 
c?r-(iV) = - f C2 - f Cl - lCs)c,C, - iCn, (197) 

/24359 172 ^ 64^2 142 n 2 1 nfAT Q^ / ,\N 



V 1620 9 5 ^ 5 



C2--Ci + ^C3)c'F\ + e{N-3){-if X 



Ci.CA|(-l)^(- 



^465 32 

- 4C2 - -cl + 54C3 - 12S__2(iV-l)C2 - 245_2,-2(iV-l) 
72 5 



" 32C3)5l(iV-l) + - 16C2)5l,l(iV-l) + y5i,i,i(iV-l) 

44 

+85i,i,2(iV-l) - —Si,2{N-l) - 8Si,2,i{N-l) + 16Si,3{N-l) 

'239 \ 88 1 1 f) 

12C2)52(iV-l) - y52,l(iV-l) + —Ss{N-l) + 8S3,i(iV-l) 

-12S4(iV-l)) + 6S-4{N-1) + 6S-4{N+1) - 12S-4{N) - —S-3{N-1) 

+^5_3(Ar+i) _ l^s_3(Ar+2) + ^S-s{N+3) + ^S^siN) - AS^3,i{N-l) 
3 5 5 3 

-45_3,i(iV+l) + 8S-3,i{N) - l5_2(iV-2) + " 4C2)5-2(iV-l) 

+ (^ - 16C2)5_2(iV+l) + f 5_2(iV+2) - - 20C2)5_2(iV) 

32 32 64 2 

+-5_2,i(A^-l) + -5_2,i(A^+l) - ^S^2.iiN) - -5„i(A^-3)C2 
3 3 3 5 

/1414 



/ ZJ. Z \ / 141:141: \ 

■(7 - 7C2)5-i(A^-2) - (-r- + 8C2 - 6C3)5_i(iV-l) 



135 
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+ 4C2 - 66C3)s_i(iV+l) - - ^C2)5_i(iV+2) - ^5_i(iV+3)C2 

+ (^ - 72C3)s_i(iV) + + 4C2)5_i,i(Ar-l) - - 28C2) 5-1,1 (iV+1) 

/724 \ 22 22 44 

+ - 32C2).S_i,i(iV) - _5_i,i,i(iV-l) - -5_i,i,i(iV+l) + y5_i,i,i(iV) 

22 22 

-45_i,i,2(7V-l) - 45_i,i,2(7V + l) + 85_i,i,2(7V) + yS_i,2(iV-l) + — 5_i,2(7V + l) 

44 

-y5_i,2(A^) + 45_i,2,i(iV-l) + 45_i,2,i(^+l) - 85_i,2,i(iV) - 45_i,3(iV-l) 

-285_i,3(Ar+l) + 325_i,3(iV) + ^Si(7V-3)C2 - ^5i(iV-2)C2 + 10Si(7V-l)C2 

78 18 4 

-105i(A^+l)C2 + -^Si{N+2)C2 - —Si{N+3)C2 - 125i(iV)C2 + -5i,_2(iV-3) 
5 5 5 

-^5i,_2(iV-2) + 20Si,_2(iV-l) - 20-Si,_2(iV+l) + ^5i,_2(iV+2) 

Of* 

-— Si,_2(iV+3) - 24Si,_2(iV) +452(7V-1)C2 + 1652(A^+1)C2 - 2052(A^)C2 
o 

+852 -2(iV-l) + 32S2 -2(A^+1) - 40^2 -2(A^)| 

{/4^7 247 ^8 8 8 

^"^^"^(^ " ^^1(^-1) - -9-^1,1(^^-1) - -5i,i,i(iV-l) + -S,,2{N-1) 

+ ^S2{N-l) + y52,l(iV-l) - y53(iV-l)) + y5_3(iV-l) + ^S.s{N+l) 

20 26 38 64 8 

--S-s{N) - -S-2{N-1) - -S-2{N+1) + -S-2{N) - -5_2,i(iV-l) 

-^5_2,i(iV+l) + y5_2,i(iV) + }^S.r{N-l) + ^5_i(iV+l) - ^5_i(iV) 

09 CO 1 nn A 

+-5_i,i(iV-l) + _S_i,i(iV+l) - — 5_i,i(iV) + -5_i,i,i(iV-l) 
+^5_i,i,i(iV+l) - ^5_i,i,i(iV) - ^5_i,2(iV-l) - ^5_i,2(7V+l) + ^5_i,2(iV)| 

+ C||(-l)^(^ + 8C2 + yd - 72C3 + 245_2(iV-l)C2 + 485-2,-2 (A^-1) 

-(y - 16C3)5i(iV-l) - (27 - 32C2)5i,i(iV-l) + 365i,i,i(7V-l) 
+485i,i,i,i(Ar-l) - 645i,i,2(Ar-l) - 365i,2(iV-l) - 485i,2,i(A^-l) + 24Si,3(iV-l) 
+ (61 - 24C2)52(iV-l) - 2452,i(iV-l) - 5652,i,i(Ar-l) + 4852,2(A^-1) 
+653(A^-1) +4853,1 (A^-1) - 1654(iV-l)) + 185_4(iV-l) + 185_4(iV+l) 



312 72 
-365_4(iV) - 32S_3(A^-1) - 765_3(A^+1) + — 5_3(A^+2) - — 5_3(iV+3) 

5 5 

+605-3(iV) - 32S-3,i(Ar-l) - 325-3,i(iV+l) + 64S-3,i(Ar) + ^S-2{N-2) 
-(^ - 8C2)5_2(iV-l) - - 32C2)5_2(iV+l) - y5_2(iV+2) 

+ (^ - 40C2)5_2(iV) + 405_2,i(Ar-l) + 565_2,i(iV+l) - 965_2,i(Ar) 
+325_2,i,i(iV-l) + 32S_2,i,i(iV+l) - 645_2,i,i(iV) - 285_2,2(iV-l) 



32 



-285_2,2(iV+l) + 565_2,2(A^) + ^5_i(iV-3)C2 + - ^C2)S-i{N-2) 
-(y - I6C2 - 12C3)5_i(iV-l) + + 8C2 - 108C3)-S_i(7V+l) 
+ (y - — C2).S_i(iV+2) + -5_i(iV+3)C2 - (101 - 96C3)5_i(iV) 

-(32 + 8C2)^-i,i(A^-l) + (84 - 56C2)S'-i,i(iV+l) - (52 - 64C2)5_i,i(Ar) 

-285_i,i,i(iV-l) - 365_i,i,i(7V+l) + 645„i,i,i(iV) - 245_i,i,i,i(iV-l) 
-245_i,i,i,i(Af+l) + 485_i,i,i,i(iV) + 325_ij,2(Af-l) + 325_i,i,2(A^+l) 
-645_i,i,2(iV) + 325_i,2(A^-l) + 325_i,2(iV+l) - 645_i,2(A^) + 245_i,2,i(iV-l) 
+245_i,2,i(iV+l) - 485_i,2,i(iV) - 205_i,3(iV-l) + 285_i,3(iV+l) - 85_i,3(A^) 

~Si{N-3)C2 + pi{N-2)C2 - 205i(Ar-l)C2 + 205i(Ar+l)C2 - ^.Si(Ar+2)C2 

Oft Q Q 

+^5i(A^+3)C2 + 2ASi{N)C2 - -5i,_2(iV-3) + -Si,.2{N-2) - 405i,_2(iV-l) 
5 5 5 

312 72 
+405i,_2(A^+l) - — Si,_2(A^+2) + — 5i,_2(A^+3) +485i,_2(A^) - 852(iV-l)C2 
o o 

-3252(iV+l)C2 + 4052(iV)C2 - 1652,_2(iV-l) - 6452,_2(iV+l) + 8052,_2(iV) 



'(A^) = (198) 



X 



Cf{Cf - ^)|-8C2 - yd - 12-S_4(iV-l) - 125_4(A^+1) + 85_3(iV-l) 

1 GO yo 

+805_3(iV+l) - — 5_3(A^+2) - — 5_3(A^+3) -405_3(Ar) + 85_3,i(iV-l) 
o o 

+8S_3,i(iV+l) + ^5_2(A^-2) - (y + 16C2)^-2(A^-1) - (y + 82(2)^-2(^^+1) 
-y5_2(iV+2) + (^ + 24C2)5-2(A^) - 85_2,i(iV-l) - 85_2,i(A^+l) 
+16S^2,i{N) + ^5_i(A^-3)C2 + (^ - ^C2)5-i(iV-2) + (^ + 20C2)5_i(iV-l) 
-(^ + 28C2).S-i(Ar+l) + (y + ^^2)s_i(iV+2) + ^5_i(7V+3)C2 
-(16 + 16C2)5_i(A^) + 165_i,i(A^-l) - 165_i,i(iV+l) - ^SiiN -3)C2 
+^SliN-2)C2 - (16C2 - 20C3)Si(iV-l) - (4OC2 - 36C3)Si(Ar+l) + ^SiiN+2)C2 
+y5i(A^+3)C2 + (32C2 - 24(3)^1 (iV) + 485i,_3(Ar-l) + 805i,_3(A^+l) 

-48Si,_3(A^) - ^5i,-2(iV-3) + ^Si,_2(iV-2) - 32Si,_2(A^-l) - 805i,_2(iV+l) 

1 68 72 

+^^i,-2(iV+2) + — 5i,_2(iV+3) + 645i,_2(iV) - 165i,_2,i(^^-l) 
5 5 

-165i,_2,i(iV+l) - 245ij(iV-l)C2 - 56^1,1 (iV+l)C2 + 485i,i(Af)C2 
-485i,i,_2(Ar-l) - 1125i,i,_2(iV+l) +965i,i,_2(iV) + 1652(iV-l)C2 
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+3252(Ar+l)C2 - 2452(iV)C2 + 3252,-2(iV-l) + 64S2,_2(A^+1) - 4852,-2(iV) 
4^'+-(iV) = 4^'+-(Ar) + 5(iV-l)| (-^ - 4C2 - f Cl + 19C3)CfCa + AC^Uf (199) 



+ 



(f + + f (I - 38C,)C^ j + i(iV-2){ + + fb)c.C, 



CfCa |45_3(7V-1) - 85_3(iV+l) + 1^5_3(Ar+2) - ^5_3(7V+3) - 165_3(Ar) 

+^5_2(iV-2) - + 4C2)5_2(iV-l) - - 12C2)5_2(iV+l) - ^5_2(iV+2) 

+ (^ - 8C2)S-2{N) + ^S_i(Ar-3)C2 + - yC2)5_i(Ar-2) 

-(^ - IOC2 - 20C3)5_i(iV-l) + + 2C2 - 60C3)5_i(iV+l) 

+(f - f C2)5_i(Ar+2) + ^5_i(Ar+3)C2 - - 40C3)5_i(iV) 

-(y - 8C2)s_i,i(iV-l) + - 24C2)5_i,i(iV+l) - - 16C2)5_i,i(A^) 

-85_i,3(iV-l) + 245_i,3(A^+l) - 165_i,3(iV) - ^5i(iV-3)C2 + Y5i(iV-2)C2 

68 18 

-145i(A^-l)C2 + 6Si{N+l)C2 - -^Si{N+2)C2 + ^5i(A^+3)C2 + 165i(iV)C2 

4 24 
--5i,_2(iV-3) + — 5i,_2(A^-2) - 285i,_2(A^-l) + 125i,_2(iV+l) 

-^.5i,_2(Ar+2) + ^5i,_2(Ar+3) +325i,_2(iV) +452(iV-l)C2 - 1252(iV+l)C2 

+852(iV)C2 + 8^2 -2(A^-1) - 24^2 -2(iV + l) + 16^2 -2(A^)| 

{8 8 16 14 62 

-5_2(A^-1) + -5_2(iV+l) - —S.2{N) - — 5_i(Ar-l) - —S.iiN+1) 

+^5_i(iV) - ^5_i,i(iV-l) - l5_i,i(iV+l) + ^5_i,i(iV)| 

{272 72 8 

245_3(iV+l) - _5_3(iV+2) + — 5_3(iV+3) + 165_3(A^) - -5_2(iV-2) 

+ (y + 8C2)5_2(7V-1) + - 24C2)5_2(Ar+l) + ^5_2(7V+2) 

/212 \ 

-(— - 16C2).S-2(A^) - 125_2,i(iV-l) - 125_2,i(A^+l) + 245_2,i(iV) 

-^5_i(7V-3)C2 - - yC2)s_i(Ar-2) + - 2OC2 - 40C3)s_i(iV-l) 

-(^ + 4C2 - 120C3)5_i(iV+l) - - l^C2)5_i(iV+2) - ^5.i(iV+3)C2 
+ (50 - 80C3)5_i(iV) + (42 - 16C2)5_i,i(Ar-l) - (78 - 48C2).S_i,i(iV+l) 
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+ ^36 - 32C2j5-i,i(iV) +85_i,i,i(iV-l) + 8S'_i,i,i(iV+l) - 165_i,i,i(7V) 

-85_i,2(iV-l) - 85_i,2(iV+l) + 165-i,2(iV) + 165_i,3(A^-l) - 485-i,3(A^+l) 
4 24 

+325_i,3(A^) + ^5i(iV-3)C2 - y5i(A^-2)C2 + 285i(A^-l)C2 - 125i(Ar+l)C2 

+— .5i(Ar+2)C2 - y.Si(Ar+3)C2 - 32Si(Ar)C2 + -5i,_2(iV-3) - — 5i,_2(7V-2) 

272 72 
+565i,_2(iV-l) - 245i,_2(Ar+l) + _5i,_2(iV+2) - y5i,-2(A^+3) - 645i,_2(A^) 

-852(iV-l)C2 + 2452(iV+l)C2 - 1652(iV)C2 - 1652,_2(iV-l) +48-S2,_2(iV+l) 
-32S2,-2(iV) 

cgJ-'^^(iV) = c(^'-"^(iV) + <^(iV-l)|(-^ + 8C2 + f C|-38C3)C^(C^-^)| (200) 
- ^) |85_3(A^-1) - 645_3(A^+1) + ^5_3(iV+2) + ^5_3(A^+3) 



+16S_3(iV) - -S_2(Ar-2) + (y + 8C2)5_2(A^-1) + + 24(2)^-2(^^+1) 

+^5_2(iV+2) - + 32C2)5-2(A^) - ^5-i(iV-3)C2 - + yC2)5-i(iV-2) 

-(^ + 20C2).S_i(iV-l) + + 28C2)5_i(iV+l) - + ^C2)5_i(iV+2) 

-^5_i(iV + 3)C2 + (16 + 16C2)5-i(iV) + ^5i(A^-3)C2 + y5i(iV-2)C2 

+ (l2C2 - 8C3).Si(Ar-l) + (36C2 - 24C3)5i(Ar+l) - — 5i(iV+2)C2 

-y5i(7V+3)C2 - (32C2 - 32C3)Si(7V) - 165i,_3(iV-l) - 48Si,_3(7V+l) 
8 32 

+645i,_3(A^) + -Si,_2(A^-3) + — 5i,_2(iV-2) + 245i,_2(iV-l) + 725i,_2(iV+l) 

128 72 
-^-5i,_2(Ar+2) - — 5i,_2(Ar+3) - 645i,_2(iV) + 165i,i(iV-l)C2 

+485i,i(iV+l)C2 - 64Si,i(iV)C2 + 325i,i,_2(iV-l) + 96Si,i,_2(iV+l) 
-128Si,i _2(iV) - 8S2(Ar-l)C2 - 24S2(Ar+l)C2 + 3252(iV)C2 - 16^2 -2(A^-1) 

-4852,-2(iV+l) + 6452,-2(Ar) 

4^'P^(iV) = <5(iV-2)|-^C^n;| + e{N-Z) (-1)^ X (201) 

f 26 
Cpufl 205_4(iV-l) + 205_4(iV+l) - 405_4(iV) + 25_3(iV-l) - —S-^{N+l) 

64 

-y5_3(iV+2) + 285_3(iV) - 165_3,i(iV-l) - 165_3,i(iV+l) + 325_3,i(iV) 

oqo 1 0S sf) 

+565_2(A^-1) - — 5_2(iV+l) + _S_2(7V+2) - —S-2{N) - 165_2,i(iV+l) 

+165_2,i(Ar+2) + 8S_2,i,i(iV-l) + 85_2,i,i(iV+l) - 16S_2,i,i(iV) - 8S_2,2(A^-1) 
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-85_2,2(iV + l) + 165_2,2(iV) + - ^C2)5_i(iV-2) - + ^C2)5_i(iV-l) 



27 ^ ' V 27 

/818 16 /448 8,\^ 88^ ,,,, 

1 04 208 208 32 

-— 5_i,i(iV-2) - _5_i,i(iV-l) + _5_i,i(iV+l) + -5_i,i(iV+2) 

1 6 28 28 

+8S_i,i(iV) + _5_i,i,i(7V-2) - _5_i,i,i(Ar-l) + _5_i,i,i(7V+l) 

1 fi 1 fi 98 98 

-^5_i,i,i(iV+2) - ^5_i,2(iV-2) + -|5_i,2(iV-l) - y5_i,2(iV+l) 

16 8 32 32 

+yS_i,2(7V+2) - -5i(Ar-2)C2 + -SiiN-l)C2 + -Si{N+l)C2 

-^5i(Ar + 2)C2 - 165i(iV)C2 - y5i,_2(iV-2) + ^5i,_2(iV-l) + y5i,_2(iV+l) 
-ySl,_2(7V + 2)-325l,_2(iV) 

4^(7V) = (202) 



'^(A^-2)|-(^-fC3)c^n^ + (^-2C3)c^n;| +e(iV-3)(-l)- x 



AT 



C^n/ 1 -105_4(iV-l) - 205_4(iV+l) + 405_4(iV+2) - 105_4(iV) + 35_3(iV-l) 

+^5_3(iV+l) - ^S^s{N+2) + ^5_3(A^+3) + yS_3(A^) + 165_3,i(iV-l) 

+ 165_3,i(^+l) - 485_3,i(A^+2) + 165_3,i(iV) + :^5_2(iV-2) 

+ 16C2)5_2(iV-l) - + 16C2)5_2(iV+l) + + 16C2).S_2(iV+2) 

"(if - 16C2)5-2(iV) - 165_2,i(iV-l) - 165_2,i(iV+l) + 725_2,i(iV+2) 
-405_2,i(iV) - 165_2,i,i(iV-l) - 85_2,i,i(A^+l) + 405_2,i,i(A^+2) - 165_2,i,i(iV) 

+125_2,2(iV-l) + 85_2,2(iV+l) - 325_2,2(iV+2) + 12S_2,2(iV) + 4'5-i(^-3)C2 

15 

+ (^ - ^C2)5_i(Ar-2) + + 24C2 + 16C3)5_i(iV-l) 

'^-^ + f C2 + 40C3)5_i(iV+l) - - ^C2 - 8C3)5_i(iV+2) - ^5_i(iV+3)C2 
+ (y - f C2 + 16C3)-S_i(iV) + (l4 + 16C2)S-i,i{N-l) - (l6 + 16C2)s_i,i(iV+l) 

-(24+ 16C2)5_i,i(Ar+2) + (26 + 16C2).S_i,i(iV) + 265_i,i,i(7V-l) 
-8S-i,i,i(Ar+l) - 725-1,1,1 (iV+2) + 545-1,1,1 (iV) + 20S-i,i,i,i(iV-l) 
-405_i,i,i,i(Ar+2) + 205-1,1,1,1 (iV) - 165-i,i,2(Ar-l) + 325_i,i,2(7V+2) 
-16S-i,i,2(iV) - 265_i,2(Ar-l) + 8S-i,2(iV+l) + 72S_i,2(7V+2) - 54S_i,2(iV) 
-24S-i,2,i(iV-l) +485_i,2,i(iV+2) - 245-1,2,1 (A^) +45-1,3 (TV- 1) 

+16S-i,3(iV+l) - 24S-i,3(Ar+2) +45-i,3(A^) - ^Si{N-S)C2 + ^Si{N-2)C2 

15 15 

-(24C2 - 8(3)^1 (iV-1) - (^C2 - 24C3)5i(iV+l) + (^C2 - 8C3).Si(iV+2) 
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-y5i(A^+3)C2 + (-^C2 - 24(3)^1 (iV) + 16Si,_3(iV-l) +485i,_3(Ar+l) 

-165i,_3(A^+2) - 485i,_3(iV) - ^Si,.2{N-3) + ^Si,^2{N-2) - 48Si,_2(iV-l) 

15 15 

fi4 Qfi Qfi 208 

- Y5i,_2(iV+l) + y5i,-2(iV+2) - —Si,-2{N + 3) + — 5i,-2(iV) 

-165i,i(iV-l)C2 -485i,i(iV+l)C2 + 165i,i(iV+2)C2 + 485i,i(iV)C2 
-325i,i,_2(iV-l) - 965i,i,_2(iV+l) + 325i,i,_2(A^+2) + 965i,i,_2(iV) 
+1652(iV-l)C2 + 1652(iV+l)C2 - 1652(A^+2)C2 - 1652(iV)C2 + 3252,-2(iV-l) 

+3252 -2(A^+1) - 32^2 -2(iV+2) - 32^2 -2(A^)| 

+ CAn/|205_4(iV-l) + 565_4(iV+l) - 7QS_4{N) + 25_3(iV-l) - ^S^3{N+1) 
388 

-— 5_3(A^+2) + 1745_3(iV) - 85_3,i(iV-l) - 485_3,i(A^+l) - WS.s,i{N+2) 

+725_3,i(Ar) + 585_2(Ar-l) - + 8C2)5_2(Ar+l) + + 8C2)5-2(A^+2) 

758 

—S-2{N) - 85_2,i(iV-l) - 45_2,i(Ar+l) + 1485_2,i(iV+2) - 1365_2,i(iV) 

+325_2,i,i(A^+l) + 165_2,i,i(^^+2) - 485_2,i,i(A^) - 325_2,2(iV+l) 
-165_2,2(iV+2) +485_2,2(iV) + " ^C2)s_i(iV-2) 

-(^ + f C2 + i4C3)5-i(Ar-i) + + f C2 + 20C3)5_,(Ar+i) 

-(^ + yC2 - 8C3)5_i(iV+2) + + I2C2 - UCs)S-i{N) - ^5_i,i(iV-2) 

-(I +4C2)5_i,i(iV-l) + + 8C2)5-i,i(Ar+l) - 1^5_i,i(7V+2) 

+ (l72 - 4C2)5_i,i(iV) + ^5_i,i,i(iV-2) - p_,,,^,{N -1) + ^5_i,i,i(Af+l) 
244 

-- 3-'5-i,i,i(iV+2) + 685_i,i,i(Ar) + 4S_i,i,i,i(iV-l) - 85_i,i,i,i(Ar+2) 

1 fi 

+45_i,i,i,i(iV) - 125_i,i,2(Ar-l) + 245_i,i,2(iV+2) - 125_i,i,2(A^) - -jS-iM^-^) 

A 28 268 

+ -5_i,2(iV-l) - y5-l,2(A^ + l) + — 5_i,2(Af + 2) - 765_i,2(iV) - 45_i,2,i(iV-l) 

+85_i,2,i(A^+2) - 45_i,2,i(iV) + 125_i,3(Ar-l) - 85_i,3(iV+l) - 165_i,3(iV+2) 
+12S_i,3(A^) - ^-Si(Ar-2)C2 + (yC2 + 2C3)Si(7V-l) - (^(2 - 12C3)-5i(iV+l) 

+ (yC2 -8(3)^1 (iV+2) - (12C2 + 6C3)5i(iV) + 85i,_3(Ar-l) + 405i,_3(iV+l) 

1 6 88 80 
-245i,_3(Ar+2) - 245i,_3(iV) - — 5i,_2(Ar-2) + _5i,_2(iV-l) - y-Si,_2(iV+l) 

+ySi,-2(iV+2) - 245i,_2(iV) - 85i,_2,i(A^-l) - 325i,_2,i(iV+l) 
+165i,_2,i(iV+2) + 245i,_2,i(iV) + 45i,i(Ar-l)C2 + 85i,i(iV+l)C2 - 125i,i(Af)C2 
+85i,i,_2(Af-l) + 165i,i,_2(A^+l) - 245i,i,_2(iV) + 852(iV+l)C2 - 852(iV+2)C2 



+1652,-2(A^+1) - 1652,-2(iV+2) 
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(2),ns 



(AT) = 



(203) 



- - yC3)c.(c. - ^) . f c|) . .(.-3) (-1)^ 

^^Fn/|-y5_2(Ar+l) + y5_2(7V) - ^5-i(iV-l) + Ms_i(Ar+l) - ^5_i(7V) 



+-5_i,i(iV+l)--5_i,i(iV) 



+ Cf(Cf - ^)|32S_3(iV+l) + ^5_3(iV+2) - ^S_3(iV+3) - 325_3(iV) 

+^5_2(iV-2) - |5_2(A^-1) - ^5_2(iV+l) - |5_2(iV+2) + ^5-2(iV) 



+|5_i(iV-3)C2 + (f - f C2)5-i(Ar-2) - - 32C2)5_i(iV-l) 



^64 32 



584 



+ (^ - - 80C3)5_i(iV+l) + (y - yC2)5_i(iV+2) + ^^5_i(iV+3)C2 



96 48 



48 



45 
'1148 



-(^ + I6C2 - 80C3)s_i(iV) + - 32C2)5_i,i(iV+l) - - 32C2)5_i,i(Ar) 

32 32 

+325_i,3(Ar+l) - 325_i,3(iV) - — 5i(Ar-3)C2 + y5i(7V-2)C2 - 325i(iV-l)C2 

-(16C2 - 16C3)5i(A^+l) - ^5i(7V+2)C2 + y5i(iV+3)C2 + (48(2 - 16C3)5i(iV) 

64 64 

+325i,_3(A^+l) - 32Si,_3(iV) - — 5i,_2(Ar-3) + — 5i,_2(iV-2) - 645i,_2(A^-l) 
96 96 

-325i,_2(iV+l) - — 5i,_2(iV+2) + — 5i,_2(iV+3) + 965i,_2(iV) - 325i,i(iV+l)C2 
o o 

+325i,i(iV)C2 - 645i,i,_2(7V+l) + 645i,i,_2(Ar) I 



184 



+ Cl\ -165_3(A^+1) + 16-S_3(Ar) - 85_2(A^-1) + — 5_2(iV+l) - — 5_2(iV) 

52 710 554 

+245_2,i(iV+l) - 245_2,i(iV) + — 5_i(iV-l) - _S_i(iV+l) + — 5_i(iV) 

+8S_i,i(iV-l) - Ms_i,i(iV+l) + yS_i,i(7V) - 165_i,i,i(A^+l) + 165_i,i,i(iV) 
+165_i,2(iV+l) - 165_i,2(iV) 



cg)^P«(iV) = <5(iV-2)|- 



80 
'27 



Cpuj \ + 0(Ar-3)(-l)^ X 



(204) 



C^n/ 1 -325_3(iV+l) + 325_3(A^) + 16S_2(A^-1) - 48S_2(-/V+1) + 32S_2(A^+2) 

+165_2,i(A^+l) - 165_2,i(A^) - ^5_i(iV-2) - ^5_i(iV-l) + ^5_i(iV+l) 

160 16 32 32 

+— 5_i(iV+2) - 16S_i(iV) - _S_i,i(iV-2) - -5_i,i(7V-l) + _S_i,i(iV+l) 

32 

-y5_i,i(A^+2) + 165_i,i(iV) 
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cg)g(iV) = 5{N-2)l -(H^ + fc3)CFnf + ^CAUf \ + 0(7V-3) (-1)^ x (205) 



f64 64 64 64 32 



-16S^2,i{N+l) + l6S.2,i{N) + ^S_i(iV-3)C2 + (§ " ^C2)5_i(A^-2) 

32 „ /32 16 ^ \ ^ /336 32 ^ \ ^ 

+ y5_i(iV-l) - (y - yC2)S_i(iV+l) + (— + yC2)5-i(iV + 2) 

-y5-i(iV+3)C2 - + Y^2)S-i{N) + 85-i,i(iV-l) - 85_i,i(iV+l) 

16 16 16 

+325_i,i(A^+2) - 325_i,i(iV) - — Si(iV-3)C2 + — 5i(A^-2)C2 + y5i(A^+l)C2 

+ y5i(Af + 2)C2 - y5i(iV + 3)C2 - Y5i(iV)C2 - ^5i,_2(iV-3) + ^5i,_2(A^-2) 

+ y5l,_2(iV+l) + y5l,_2(A^ + 2) - y5l,_2(iV + 3) - ^5l,_2(iV)| 

CAn/|-965_3(A^+l) + 965_3(A^) + 165_2(A^-1) - 805_2(A^+1) + 2085_2(iV+2) 

-1445_2(iV) + 645_2,i(iV+l) + 325_2,i(iV+2) - 965_2,i(iV) - — 5_i(iV-2) 

9 



32 32 

-—S.iiN-1) + — S_i(iV+l) - + 16C2j5-i(iV+2) + (96 + 16C2j5-i(iV) 

-^5_i,i(iV-2) - ^5_i,i(iV-l) + ^5_i_i(Ar+i) _ ,(iv+2) 
+1605_i,i(iV) - 325__i,i,i(7V+2) + 325_i,i,i(iV) + 325_i,2(A^+2) - 325_i,2(A^) 
-325i(Af+l)C2 + 16Si(iV+2)C2 + 16Si(iV)C2 - 645i,_2(iV+l) + 325i,_2(A^+2) 

+325i,_2(iV) 

Appendix C 

Here we present the x-space expressions of the 2-loop coefficient functions. We find 



3 



cg'+-(x) = (206) 

^ ^ f 3229 4 2191 36 o , , x/ 3155 11 1 , , 

CfCa< \ X + 4pgJx){ \ C2 + -Cs - 3H_2oix) 

\ 180 5x 20 5 432 3 2^ ^ 

-^i?o(x) + Fo(x)C2 - ^Ho,o{x) - ^Fo,o,o(x) - ^^i(^) + 3Fi(x)C2 

-iii7i,o(x) - 2Fi,o,o(x) - ^^i-i(^) + ^1,1,0(2;) - ^1,2(2;) - ^H2{x) - Hsix] 
6 6 3 

-^(133 - 3nx)Hi{x) - 4(5 + + 2; - 6^2 - ^x3)f_i,o(x) + :^(13 + — 
\ 5x^ 5 ^ 30 x 

+17532; - 216x^)Ho{x) + 4(1 - 5x)(C3 + Hi{x)C2 - H_2,o{^) - Hi,o,o{x 



3 



The expression for the coefficient function ^ in eq. (13) of ref.|l|] contains a typographical error. The term 



+488/27a;2 should read +448/27a;^. 
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+2(l -3x + 12a;2 + y^^)'^^ - 4(i - x + Gx^ + -x^^Ho,o{x) - 4(1 + x)H2{a 

il-x) 



5 ^' V 5 

/5465 251 71 2 140. 

{ 23 27 4 /247 19 5 29 

+ 2iJ2(a;)) - ^(7 + 19x)/fo(x) - ^(1 + lSx)Hi{x) 

/457 38^ , .1 

+6/fo(x)C2 - ^^0,0(2;) - 2/^0,0,0(2;) + yi?i(2;) + 2Hi{x)C2 - Ihi,o{x) 
-3i/i,o,o(x) - ^Hi^iix) - 8i/i,i,o(a;) - 6i7i,i,i(x) - 6/7i,2(x) - 3if2(x) - 6H2,o{x) 

+8(5 + -—^ + x-6x^- -x^)i/_i,o(a;) + (29 + 25x + 48^^ + —x^^Hoflix) 
-12(1 - 3x)C3 + 8(1 - 5x)(ii-_2,o(x) + i?i,o,o(^c) - Hi{x)C2) 
+2(1 + x)(5Ho,o,o{^) - 4^o(x)C2 + 7Hi,o{x) + 2H2,o{x) + 2H2,i{x) + AH^{x)) 
-24(1 + X + 2x2 + -x3)C2 + 4(7 + llx)if2(x) + 2(5 + 9x)i7i,i(x) 

+ (^ + 69C2 + 6C| - 78C3)<5(1 - x)| , 

c?}f-\x) = (207) 

Cf{Cf - ^) + a + + + 4pqq(-x) (7(3 + 6H^2.o{^) ~ 8if_i(x)C2 

-8i7_i,_i,o(a;) + 10i/_i,o,o(a;) + 4i/_i,2(x) - 2R^{x) + 2i/o(a;)C2 - 3ffo,o,o(a;) 
-2i73(x)) + 32(1 + -1^ + a; + - A^3jiy_^ o(^) _ 16(1 - ar)/fi(x) 

--(13 + - + 53x - 36x2 Wo(x) _ g/^ ^ x)H2(x) - 8fl + 4x + Gx^ - -x^Woo(a;) 
5 X / V 5 ^ ' 

+8(1 + 5x)(C3 + H_2,o{x) - 2i/_i,_i,o(x) + i?-i,o,o(x) - i^-i(a;)C2) 
+4(1 + 7x + 12x2 - Y^^)C2| , 

+ ^ - + yx2 + ^(25 - 71x)ifi(x) - 8(1 - 3x)(C3 + H,{x)C2 

-i?_2,o(x) - i?i,o,o(a;)) + 4(6 + ^ + ^ + 2x - 5x2 - ^^^^^H.^^oix) 
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-^(31 + ^ + 121a; - 27x^)Ho{x) + a(i + 3x - hx"^ - ^x^) (C2 - Hq,q{x)) | 
+ CFnf\^-^ + ^x + ^{l + x){2Ho{x)+H,{x))^ 

+ C J-^ - A + 1^^ _ 11^2 _ 2(21 _ 39x)ili(x) + 1(2 + - + %- 9x')Ho{x) 
[5 ox 55 5^x2 / 

-Sfe + + - + 2x - 5x2 - 9 3X ^ ^/^ _ ^ ^^^2 ^ ^a;3\^2 

V X 5 / V 5 ^ 

+16(1 - 3x)(C3 + Hi{x)C2 - H_2,o{x) - Hi,o,o{x)) 

-8(1 + x)(Hi,o{x) + Fi,i(x) + Ih2{x)) + 8(2x - 5x2 _ ^^s^Hoflix)^ , 

= 4^-^^(x)+ (209) 
/ Ca\ f 312 8 232 72 9 / 11 9 9 

+ ^ (7 + ^ + 37x - 9x2] ^^(^^ + 8(^1 _ _ 5^2 ^ 9^3^ j^^^ _ q(^)] 
-16(1 + 3x)(C3 - /^-i(x)C2 + H^2,o{x) - 2iJ_i,_i,o(x) + iJ-i,o,o(x)) I , 

4^'P^(x) = (210) 

^ fl58 344 422 448 o 8 / ^ 13 4 2\ ,w x 

"^^n^ + 27^ - + + 3 - 3^ - + r 

-16(l + ^ + a; + Ix^)H-i,o{x) + 4(1 + ^ - X - ^x2) (i/i,o(x) + H^,i{x 
-16x^H2{x) - 2(1 - 15x + ^x^)Ho^oix) + 56(l - ^x - ^x^)Ho{x) 
-8(1 + x)(C3 + 2/fo(^)C2 - |i?OAo(a^) - H2,o{x) - H2,i{x) - 2H3{x)) 



-f(i + 3x-3x2)c.), 



cg(x) = (211) 



^ f 647 8 239 36 2 1 



4 2 

nfCF{ + — + ^x - ^x^ + 48( 1 + + 9^ + ;^a;^)^^-i,o(x) 



-4(59 + - - ^x + 162x2)i7o(x) - 3(1 - ^x + 24x2 + '^Ix'^R^fsix) 
15^x4 / v9 5^ 



-8(2 - 7x + 9x2)i/2(x) - 2(13 - 40x + 36x2) (^^^ ^(3.) + ^^^^(^ 

13 9 6 

-2(7 - 20x + 12x2)i/i(x) + 16(l --^x + -x2 + -x3)C2 + 8(4 + 9x2)C3 
+32(1 + x2)i7_2,o(x) + 16x2(i7_i(x)C2 + 2i7_i,_i,o(x) - H-x^^M + ^o(^)C2 
-^/fo,o,o(^) + iix{x)C,2 - Hi,oAx) - ^H2,o{x) - ^H2,i{x) - H3{x)) 
-16(iJ_i(x)C2 + 2iJ_i,_i,o(x) - H_i,o,oix))p^g{-x) + 16(Ho{x)C2 - ^^^OAo(aj) 
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11 5 3 3 



^ [239 344 1072 4493 o 4 67 2\ / ^ 

+ n,C.|- + ^ + - - 4(1 - ^ - 20x + -x^)H,,oix) 

-4(1 - ^ - 18. + |x>i,i(x) + 8(1 - ^ - 18x + 



-2(1 - 88a; + —x'^^Ho^oix) - 8(l - 18x + —x'^)H2{x) + a(i - 12x + Qx^)c,z 

/ 292 1045 oN , , 
+4(5 + Ux)HoM^) + 58(1 + ^x - -^x^)Ho{x) 

-8(1 + 8a; - 2x^) [Ho{x)C2 - H^ix)) + 16a;(3 - x)(^H2,o{x) + H2,i{x)) 
+lQx'(H_2fi{^) - ^i?-i(a;)C2 - H^r,-ifl{x) + ^if-i,o,o(^) - \hi{x)C2 

+\Hi^Ofi{x)) - 4(if_i(x)C2 - 2i/_i,_i,o(x) - 2i7_i,o,o(x) - 2i/_i,2(x))pqg(-x) 

3 3 11 1 

+8(/fi(x)C2 - 2^i,o,o(x) - 2^i,i,o(^) - 2^i.i,i(^) - 2^i.2(a;))pqg(a;) I , 



„(2),ns. 



^lT = (212) 
""^^^{^ " " ^^(2^o(x) + i?i(a;)) I 

^ Cyi\f392 64 6916 96 n / 2 3 ox , , 

+ C.(C. - -) |- + - - + -x^ + 32(2 + ^ + X - -x3)i/_,o(x) 

+f (1 - ^ - f X + 3x2)iJo(a;) - ^xi7i(x) + 32x(l - ^-x") (c,2 - i?o,o(x)) 
+32x(2C3 + C,2Hi{x) - C2H-i{x) - 2iI_i,_i,o(x) + i?-i,o,o(a^) - H^,Ofi{^)) | 

+ C-^j-f + - 8(1 - f ^)^o(x) - 8(1 - ^x)i^i(x) - 16x(^C2 - Ho,o{x) 



-Hi,o{x)-H,,,{x)-^H2{x 
=S''(x) = (213) 

-16x(C2-2/fo,o(a;)-/f2(a;))|, 

4'i(^) = (214) 

^ f 128 32 304 336 2 8 / 4 „ 2\ . x 

"/C^F*^ -- + + - — (13 + - + 78x - 36x^)Ho{x) 

' y 15 15x 5 5 15 V x J ^ ' 
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-8(l + 3x- Ax^)Hi{x) + Y^(^ - + Qx^)h_i,o{x) - yx(l + ^x^)Ho,o{ 



16 12 2 



+Yx(^1 + y2;^jC2 - 16xi72(x)| 

+ ^/C'aI^ - ^ + ^x - ^x2 + (3 - i + 27x - 29x2)i/i(x) - 32x(2 - x)C2 
13 9x 3 9 3^x ^ 

+16 (l + 8x - 13x2)i?o(x) + 32x(l - x)(iIi,o(x) + i?i,i(x)) + 32x(3 - x)i?2(x) 
+32x(l +x)i/_i,o(x) + 96xi/o,o(a;)| . 
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